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1. 


SUMMARY 


Recent  work  on  earth  penetrating  devices has  focused  on  the 
terradynamlcs  of  the  target  material  and  on  the  structural  response  of 
the  penetrator.  Test  results^"^  show  that  laige  bending  strains  can 
be  produced  in  the  penetrator  casing  when  i 3 Impacted  at  an  angle  of 
attack  of  only  a few  degrees.  V)hen  these  stialns  exceed  the  yield 
strain,  local  weakening  of  the  cross  section  leads  to  breakup  of  the 
penetrator;  thus  it  Is  desirable  to  keep  the  penetrator  response  elastic. 
Because  this  response  la  sensitive  to  both  Impact  velocity  V and  angle 
of  attack  a,  a tradeoff  in  these  two  Impact  parameters  can  be  made  In 
the  design  of  an  earth  penetrator.  The  primary  objective  of  our  work 
was  to  determine  this  tradeoff  In  the  form  of  critical  Impact  curves  for 
representative  penetrator  structures.  A second  objective  was  to  demon- 
strate that  an  explosive  loading  technique  can  be  used  to  simulate 
angle-of-attack  Impacts  in  the  laboratory. 

Construction  of  critical  Impact  curves  requires  knowledge  of  both 
the  structural  response  of  the  penetrator  and  the  response  of  the  target 
material  (terradynamlc  response) . In  this  work  these  two  responses  were 
treated  Independently,  then  combined  by  matching  the  terradynamlc  loads 
with  the  loads  applied  to  the  penetrator  structure. 

Structural  response  was  calculated  with  a mathematical  model  based 
on  elastic  Timoshenko  beam  theory  with  axial  thrust.  This  model  allows 
efficient  calculation  of  the  elastic  response  of  penetrators  under  com- 
bined axial  and  lateral  loading.  The  peak  response  stress  was  calculated 
for  a range  of  load  rise  times  and  axial  and  lateral  load  amplitudes. 

Terradynamlc  response  and  the  associated  loads  that  are  applied  to 
the  penetrator  structure  are  not  easily  determined.  However,  available 
experimental  data  and  computer  code  calculations  show  that  the  loads 
have  three  characterlotlcs . First,  the  shape  of  the  resultant  load 


history  rises  llneraly  to  a peak,  followed  by  a constant  value  where 
the  rise  time  Is  the  nose  length  divided  by  the  Impact  velocity.  Second, 
the  peak  loading  force  Is  proportional  to  impact  velocity.  Third,  based 
on  comparison  of  calculated  structural  response  and  of  response  measured 
In  an  angle-of-attack  reverse  ballistics  test,  the  angle  between  the 
resultant  force  and  the  penetrator  axis  Is  three  times  the  angle  of 
attack. 

With  these  three  load  characteristics,  the  peak  stress  calculated 
as  a function  of  load  rise  time  and  amplitudes  can  be  expressed  In 
terms  of  the  Impact  velocity  and  angle  of  attack.  We  then  can  construct 
critical  Impact  curves  that  show  the  combinations  of  V and  a that 
produce  a given  level  of  peak  strain,  or  factor  of  safety,  In  the  structure. 
The  prlnclpa]  advantage  of  critical  Impact  curves  Is  that  they  characterize 
the  response  over  a range  of  Impact  conditions.  For  example,  critical 
impact  curves  for  four  structures  are  shown  In  Figure  1.*  These  curves 
can  be  used  to  select  the  penetrator  material,  on  the  basis  of  yield  stress, 
for  a given  range  of  V and  a within  which  the  penetrator  must  function. 
They  can  also  be  used  to  interpret  results  of  experiments  or  Interpret  more 
detailed  load  and  structural  response  calculations  (e.g.,  a finite  element 
code  prediction).  For  example,  the  appropriate  values  of  V and  a can 
be  selected  to  minimize  the  tests  or  code  calculations  needed  to  define 
the  curve. 

The  critical  Impact  curves  are  perhaps  most  useful  for  making  design 
tradeoffs  among  candidate  penetrator  structures.  For  example,  the  curves 
for  Structure  B (deep  penetrator)  and  for  Structure  C (shallow  penetrator) 
pass  through  a similar  region  In  the  V-a  plane.  Both  these  penetrators 
can  be  made  stronger  by  Increasing  only  the  wall  thickness  (from  that  of 
Structure  C to  that  of  Structure  A)  or  by  decreasing  the  length  (from 
that  of  Structure  B to  that  of  Structure  A).  For  small  angles  of  attack, 


^1000  average  pressure  over  the  frontal  area  of  the  penetrator 

at  an  Impact  velocity  of  1000  ft/sec  and  is  a characteristic  of  the 
target  material. 


10 


however,  the  velocity  range  can  be  extended  more  by  Increasing  wall 
thickness  than  by  decreasing  length.  Similar  conclusions  can  be  drawn  by 
comparing  curves  for  Structures  B and  C with  that  of  Structure  D. 

To  meet  the  second  objective  of  this  work,  demonstration  of  a 
technique  for  simulating  angle-of-attack  Impacts,  we  used  an  existing 
explosive  simulation  device.  This  device  produces  resultant  axial  and 
lateral  time-varying  loads  similar  to  those  occurring  In  angle-of-attack 
Impacts,  but  not  the  detailed  load  distribution.  The  controlled  flow  of 
high  pressure  gases  from  a confined  explosion  is  used  to  load  a piston 
that,  in  turn,  loads  a penetrator  initially  at  rest.  Since  we  simulate 
only  the  damage-producing  portion  of  the  load,  which  extends  up  to  and 
slightly  beyond  the  maximum  load,  the  kinetic  energy  Impacted  to  the 
penetrator  Is  much  less  than  that  required  in  ballistics  tests;  thus, 
simulator  tests  can  be  performed  conveniently  in  the  laboratory.  This 
technique  also  allows  hard-wired  measurement  of  the  structural  response 
of  the  penetrator. 

Several  penetrator  structures  were  tested  in  nominally  1/4  scale. 

For  example,  the  model  typical  of  deep  earth  penetrator  structures  (thick- 
walled  model)  is  made  of  AISI  1020  steel,  Is  9 Inches  (22.9  cm)  long,  and 
has  a 0 .87S-lnch-dlameter  (2.22-cm-dlameter)  cylindrical  cavity  over  the 
aft  two-thirds  of  its  length.  The  load  produced  on  this  penetrator  at 
an  angle  of  attack  was  simulated  by  tilting  the  model  penetrator  through 
an  angle  of  19“30'  with  respect  to  the  penetrator  axis.  The  peak  axial 
compressive  strain  was  0.022  percent.  The  peak  total  (axial  plus  bending) 
compressive  strain  was  about  0.042  percent.  This  bending  response  is 
similar  to  that  observed  in  angle'-of-attack  reverse  ballistics  tests 
performed  by  AVCO.  '’ 

These  and  other  test  results  show  that  the  load  simulator  can  apply 
to  model  penetrators  loads  similar  to  those  occurring  in  angle-of-attack 
impacts.  Ttie  tests  also  indicate  that  the  loader  could  be  built  in  a 
larger  size  to  test  full-scale  penetrators. 


2. 


INTRODUCTION  AND  DETAILED  SUMMARY 


BACKGROUND 

Earth  penetrating  devices  have  potential  for  use  as  tactical 
weapons.  Such  a device  would  have  a projectile  shape  and  sufficient 

•N 

impact  velocity  to  enter  soil,  soft  rock,  or  concrete  and  travel  into 
the  target  before  a charge  in  the  ponetrator  is  detonated. 

A variety  of  penetrator  structures  have  been  proposed  for  different 
impact  conditions.  For  deep  penetration  into  hard  targets,  penetrators 
with  a radluB-to-thickness  ratio  of  about  3 and  a length-to-dlameter 
ratio  of  about  10  have  been  used.^  A thick  wall  was  chosen  to  withstand 
the  impact  stress,  and  a long  slender  shape  was  chosen  to  minimize 
broaching.  For  shallow  penetration  (as  little  as  one  penetrator  length) 
into  soft  targets,  a radius-to-thickness  ratio  of  about  10  and  a 
length-to-dlameter  ratio  of  3 to  6 have  been  proposed. 

Determining  the  feasibility  of  deploying  an  earth  penetrating 
weapon  requires  investigation  into  three  areas  associated  with  impact 
and  penetration.  First,  the  early-time  Impact  response  of  the  casing 
must  be  understood  so  that  it  can  be  designed  to  stay  intact  and  allow 
penetration  of  the  target.  Second,  the  response  of  the  internal 
components  must  be  understood  so  that  they  can  be  made  to  function  after 
the  impact.  Third,  the  terradynaraics , or  motion  of  the  penetrator 
through  the  target  media,  must  be  understood  so  that  the  system  can  be 
designed  to  penetrate  to  the  required  depth  for  detonation. 

Most  of  the  previous  work  on  earth  penetrators  has  focused  on 
terradynamics. These  investigations  have  demonstrated  experi- 
mentally that  solid  or  thick-walled  projectiles  can  travel  tens  of 
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feet  In  soil  nnd  soft  rock.  In  addition  they  have  yielded  analytical 
techniques  for  calculating  the  forces  on  a rigid  penetrator  and  the 
resulting  motion  of  the  penetrator  and  the  target  for  a rigid  pene- 
trator under  normal  Impact. 

Both  calculated  Impact  loads^  and  measured  acceleration  response^ 
show  that  the  resultant  force-history  for  normal  Impact  consists  of  two 
distinct  parts t as  Figure  2 Illustrates:  an  approximately  linear  rise 
to  a peak  over  the  time  t^  required  for  the  structure  to  penetrate  to 
Its  full  diameter  I and  a very  gradual  decay  associated  with  the  rigid 
body  deceleration  of  the  structure.  For  a nominal  full-scale  structure 
[6  Inches  (15.24  cm)  In  diameter,  60  Inches  (152.4  cm)  long,  and 
weighing  400  pounds  (181  kg)]  Impacting  sandstone  at  1500  ft/sec 
(457  m/sec),  the  loading  rise  time  Is  about  1 msec  and  the  peak  force 
Is  about  500,000  pounds  (224  kN).‘^ 

Less  work  has  been  done  on  determining  the  response  of  penetrator 
structures  or  Internal  equipment.  In  particular,  we  need  to  know  the 
loads  that  produce  damage  In  a penetrator  casing  and  how  these  loads 
depend  on  the  structural  or  loading  parameters  that  are  at  the  disposal 
of  a designer.  Nevertheless,  existing  results  from  tests  on  deep- 
penetrator  structures  have  Identified  some  of  the  important  parameters. 

In  reverse  ballistic  tests  performed  by  AVCO,^  relatively  low  strains 
were  produced  In  a simple  penetrator  structure  under  a normal  impact 
with  a rock  simulant;  however,  under  an  angle-of-attack  Impact,  much 
larger  strains  were  produced  by  the  bending  Induced  In  the  penetrator. 

In  ballistic  tests  performed  on  scale  model  penetrators  by  Martin 
Marietta  Aerospace  Company,^  the  penetrator  structure  failed  in 
angle-of-attack  Impacts;  these  tests  results  also  Indicate  that  large 
strains  and  failure  were  caused  by  bending. 

These  observed  failures  In  penetrators  under  angle-of-attack  Impacts 
are  postulated  to  occur  as  follows.  The  axial  component  of  the  load 
produces  a compressive  stress  along  the  entire  length  of  the  penetrator. 
The  lateral  component  of  the  load  produces  bending  stresses  whose 
magnitudes  In  tension  and  compression  are  equal  at  a given  axial 
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location.  Thus  the  total  stress  of  greatest  amplitude  Is  compressive 
and  produces  local  yielding  of  the  penetrator  wall.  The  weakened 
cross  section  at  this  location  then  allows  continued  deformation  and 
ultimately  fracture  under  the  large  tensile  strains  at  this  location. 

Thus,  Incipient  yielding  In  compression  defines  the  onset  of  a divergent 
Instability.  The  response  that  leads  to  yielding,  that  Is,  the  elastic 
response,  therefore  determines  If  failure  occurs.  Thus  the  analysis 
need  predict  only  the  elastic  response  and  the  critical  response  para- 
meter Is  the  maximum  total  compressive  stress. 

Also,  observed  maximum  compressive  strains  and  failures  In  pene- 
trator structures  occur  at  distances  greater  than  one  diameter  from  the 
penetrator  nose.®”®  The  stress  distribution  over  the  cross  section  at 
these  locations  depends  only  on  the  resultant  forces  applied  to  the  end 
of  the  structure.  Thus  failure  depends  only  on  the  resultant  axial  and 
lateral  loading  forces  with  the  appropriate  time-history,  and  not  on  the 
details  of  the  pressure  load  distribution  on  the  penetrator  nose.  We 
will  see  that  this  observation  allows  us  to  use  a relatively  simple 
experimental  technique  and  response  analysis. 

OBJECTIVES  AND  APPROACH 

The  primary  objectives  of  this  work  were  to  develop  a simple  analysis 
that  accounts  for  penetrator  bending  induced  by  angle-of-at tack  impacts 
and  to  determine  the  effect  of  structural  and  loading  parameters  on  the 
survivability  of  a penetrator  casing  on  impact.  These  results  can  then 
be  used  to  (1)  Identify  those  parameters  to  which  the  damage  producing 
loads  are  sensitive,  e.g.,  impact  velocity  and  angle  of  attack; 

(2)  provide  a means  of  arriving  at  an  optimum  structural  design  within 
the  par.^imeter  ranges  dictated  by  system  requirements;  and  (3)  select 
specific  configurations  and  loads  for  testing  and  for  more  detailed 
calculations . 


A secondary  objective  was  to  demonstrate  the  capability  of  an 
explosive  load  simulator,  already  developed  on  a separate  DNA  contract,* 
for  loading  a scale  model  penetrator-llke  structure  with  a resultant 
load  simulating  an  Impact  load.  These  experimental  results  were  used 
to  validate  the  analysis  developed  here.  Also,  by  developing  and  demon- 
strating the  usefulness  of  this  device  in  small  scale,  we  have  made  the 
concept  readily  available  for  testing  full-scale  penetrator  structures, 

The  analytical  approach  to  predicting  penetrator  response  Is  based 
on  elastic  Timoshenko  beam  theory  with  axial  thrust.  This  engineering 
theory  treats  the  loading  and  response  In  terms  of  the  resultant 
longitudinal  force,  transverse  force,  and  bending  moment  at  any  cross 
section  along  the  structure.  The  nose,  aft  mass,  and  payload  are 
modeled  as  rigid  masses.  A one-dimensional  finite  difference  character- 
istic grid  along  the  length  of  the  structure  can  then  be  used  in  the 
numerical  solution  of  the  governing  equations.  This  procedure  is  an 
efficient  way  to  calculate  the  elastic  response  of  penetrators,  allowing 
a broad  range  of  loading  and  structural  parameters  to  be  examined  at  a 
low  cost.  It  is  not  intended  to  replace  more  elaborate  finite  element 
and  finite  difference  codes,  but  rather  to  supplement  code  calculations 
in  the  design  stage  when  a number  of  configurations  are  being  considered 
and  structural  details  are  yet  to  be  determined. 

The  impact  loads  developed  and  their  dependence  on  the  impact 
conditions  are  not  well  understood,  especially  under  angle-of-attack 
impacts.  Therefore,  the  load  is  treated  as  a parameter;  that  is,  we 
calculate  the  response  for  a range  of  loading  parameters  as  well  as  for 
the  parameters  that  describe  the  penetrator  casing.  The  loads  are  re- 
lated to  the  impact  conditions  through  available  experimental  and 
analytical  load  data. 
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The  experimental  approach  makes  use  of  an  existing  explosive 
simulation  device  that  produces  resultant  axial  and  lateral  time-varying 
loads  similar  to  those  occurring  In  angle-of-attack  impacts,  but  not  the 
detp'*Led  load  distribution.  The  device  uses  the  controlled  flow  of 
high  pressure  gases  from  a confined  explosion  to  produce  loads  on  a 
penetrator  initially  at  rest.  This  technique  allows  hard-wired  measure- 
ment of  the  structural  response  of  the  penetrator  and  avoids  the  high 
velocities  required  to  produce  loads  by  target  impact  (we  simulate  only 
the  damage-producing  portion  of  the  load  that  extends  up  to  and  slightly 
beyond  the  maximum  load) . 

The  experiments  are  summarized  first.  Then  we  describe  the  develop- 
ment of  the  analysis,  comparison  of  predicted  and  measured  responses, 
and  the  application  of  the  analysis  to  determining  critical  Impact 
curves . 

EXPERIMENTS 

Figure  3 Is  a sectioned  assembly  drawing  of  the  loading  fixture. 

As  shown,  It  Is  a configuration  for  simulating  normal  impact  loads  on 
1/4-scale  model  penetrators.  High-pressure  gaseous  explosive  products 
flow  through  the  orifices  and  transmit  the  load  to  the  penetrator 
through  the  piston.  The  rise  time,  duration,  and  decay  time  can  be 
varied  by  using  different  size  spaclngs  and  vent  holes. 

To  simulate  the  load  on  a penetrator  that  Impacts  at  an  angle  of 
attack,  the  piston  Is  designed  to  produce  a combined  axial  and  lateral 
loading  and  to  measure  the  load  applied  to  the  penetrator.  The  com- 
bined loading  la  produced  by  tilting  the  penetrator  through  an  angle  6 
with  respect  to  the  piston  axis.  During  the  loading  the  penetrator  is 
allowed  to  slide  relative  to  the  piston.  The  resultant  vertical  force 
history  Is  measured  with  a load  cell  Inside  the  piston;  the  resultant 
horizontal  force  history  is  deduced  from  the  measured  motion  of  the 
penetrator. 
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FIGURE  3 ASSEMBLY  DRAWING  OF  TEST  FIXTURE  FOR  SIMULATING  IMPACT  LOADS 


All  the  model  structures  tested  have  a 1.50-inch  (3.81-cm)  outer 
diameter  and  a hemispherical  front  end,  rather  than  a pointed  nose, 
to  allow  the  desired  load  to  be  applied  more  easily  and  accurately. 

The  model  typical  of  deep  earth  penetrator  structures  (thick-walled 
model)  is  made  of  AISl  1020  steel.  Is  9 inches  (22.9  cm)  long,  and  has 
a 0.873-lnch-dlameter  (2. 22-cm-diameter)  cylindrical  cavity  over  the 
aft  two-thirds  of  its  length.  Each  model  is  instrumented  with  eight 
axial  strain  gages,  four  at  the  front  station  about  2 inches  from  the 
nose  and  four  at  the  aft  station  at  about  the  midlength.  At  each 
station  the  gages  are  uniformly  spaced  around  the  circumference  and 
oriented  so  that  the  gages  measure  only  strain  in  the  axial  direction. 

Figure  4 shows  the  load  cell  and  strain  records  from  Test  64  of 
the  thick-walled  model.  The  initial  angle  of  tilt  6 was  19*30* 

(0.340  rad).  In  this  test  the  peak  load.  Figure  3(a),  was  13,000  pounds 
(66.7  kN) . The  peak  axial  compressive  strain.  Gages  5 and  7,  Figure  3(b), 
was  0,022  percent.  The  bending  strain  reduced  the  strain  at  Gage  6 but 
added  to  the  strain  at  Gage  8 to  produce  a peak  compressive  strain  of 
about  0.042  percent.  This  bending  response  is  similar  to  that  observed 
in  angle-of-attack  reverse  ballistics  tests  performed  by  AVCO.®  At 
higher  load  levels,  yielding  would  first  occur  in  the  vicinity  of  Gage  8, 
followed  by  a reduction  in  bending  stiffness,  buckling,  and  fracture  of 
the  penetrator  structure. 

These  and  other  test  results  show  that  the  load  simulator  can  apply 
to  model  penetrators  loads  similar  to  those  occurring  in  angle-of-attack 
Impacts.  The  structural  response  of  the  penetrator  has  been  accurately 
measured,  using  hard-wired  strain  gages.  The  tests  indicate  that  the 
loader  could  also  be  built  in  a larger  size  to  test  full-scale 
penetrators. 

ANALYSIS  OF  PENETRATOR  RESPONSE 

In  the  analysis,  the  central  portion  of  the  penetrator  Is  modeled 
by  elastic  Timoshenko  beam  theory  including  thrust.  This  formulation 
Includes  the  two  dominant  response  mechanisms  governing  angle-of-attack 
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FIGURE  4 TEST  64-19"30'  (0,340  rad)  ANGULAR  LOADING  OF 
THICK-WALLED  MODEL 
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Impacts:  axial  compression  and  bending  under  lateral  loads.  The  front 
portion  of  the  solid  nose  and  the  solid  aft  end  are  modeled  as  rigid 
masses.  Loading  of  this  structure  Is  specified  by  the  magnitude,  time 
history,  and  orientation  of  a resultant  Impact  force  through  a fixed 
point  in  the  front  mass. 

The  response  predicted  by  the  analysis  was  first  compared  with  the 
penetrator  strains  measured  In  the  simulator  experiments.  The  measured 
loads  were  used  as  input  for  the  analyses.  To  account  for  the  small 
uncertainty  in  the  point  of  application  of  the  resultant  loading  force 
in  the  simulator  experiments  (due  to  deformation  of  the  face  of  the 
loading  piston),  upper  and  lower  bounds  on  the  response  were  calculated 
for  extreme  locations  of  the  point  of  load  application  for  each  experi- 
ment. Thus,  these  comparisons  not  only  provided  a check  on  the  analysis, 
but  also  aided  in  the  understanding  of  the  simulator  loading  technique. 

For  the  thick-walled  penetrator  of  Teat  64,  Figure  5 shows  the 
upper  and  lower  bounds  on  calculated  strains  along  with  the  strain 
measured  at  strain  Gage  8.  Hie  error  bands  on  each  of  the  predicted 
strains  correspond  to  the  uncertainty  In  the  horizontal  force  and  in  the 
angle  of  inclination.  The  total  strain  calculated  for  the  two  extreme 
points  of  load  application  bound  the  initial  peak  measured  in  the 
experiment . 

The  analysis  was  then  used  to  investigate  the  effects  of  the  ratio 
of  lateral  to  axial  load,  pulse  shape,  rise  time,  end  masses,  and  lateral 
payload  inertia.  We  discuss  here  only  the  effect  of  rise  time.  This 
elfect  was  studied  by  applying  loads  with  the  different  rise  times  shown 
in  Figure  6(a)  and  with  a ratio  of  lateral  to  axial  load  of  0.2.  As 
discussed  in  Section  5,  this  corresponds  to  an  angle  of  attack  of  about 
3.8  degrees.  The  structure  was  a simple  steel  tube  with  length-to-dlameter 
ratio  ltd  ■ 6 and  radius-to-thickness  ratio  a/h  ■ 4. 

The  bending  and  axial  stress  histories  at  Station  x ■ x/f.  » 0.4828, 

normalized  with  respect  to  the  normal  stress  o at  the  loaded  end  of 

8 

the  tube  are  shown  in  Figure  6(b). 
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The  largest  bending  stress  * ±1.51)  Is  produced  by  the 

loading  of  zero  rise  time  (step  Input).  In  order  of  Increasing  loading 
rise  times,  the  peak  bending  stresses  a. /a  produced  are  1.48,  1.40, 
1.31,  and  1.20.  As  the  rise  time  Increases  further,  the  peak  stress 
at  this  location  approaches  the  steady-state  value  of  **  0.793. 

Furthermore,  the  time  at  which  the  peak  occurs  Increases  with  rise  time. 

For  the  axial  stress,  the  peak  stress  at  a given  location  Is 
determined  by  whether  the  loading  rise  time  Is  less  than  or  greater 
than  the  time  at  which  a reflected  tensile  stress  wave  from  the  aft 
free  end  arrives  at  the  location.  For  loading  curves  with  rise  times 
of  leas  than  one  transit  time,  the  maximum  axial  stress  produced  Is 
a /a  ■ -1.00.  For  the  loading  curves  with  rise  times  of  more  than  one 

Bi  S 

transit  time,  the  maximum  stress  Is  a /a  > -0.617.  The  steady-state 

a s 

value  at  this  location  Is  a /a  " -0.517. 

a s 

Thus,  wo  conclude  that  the  largest  effects  due  to  Increasing  the 
loading  rise  time  are  the  Increase  In  the  times  at  which  peak  bending 
occurs  and  the  variation  In  axial  stress  for  loading  rise  times  near 
one  transit  time.  The  variation  In  loading  rise  time  has  less  effect 
on  the  peak  bending  stress  and  on  the  peak  axial  stress  for  rise  times 
larger  than  two  transit  times. 

CRITICAL  IMPACT  CURVES 

In  this  section  we  apply  the  analysis  to  some  of  the  penetrator 
struitures  of  Interest.  These  structures  range  from  deep  penetrators 
(large  i/d,  small  a/h)  to  shallow  penetrators  (small  Z/d,  large  a/h). 
Calculations  were  also  made  for  penetrators  with  Intermediate  values  of 
Z/d  and  a/h  so  that  we  could  see  how  these  ratios  affect  response. 

For  design  purposes  the  most  useful  Information  is  the  relationship 
between  the  Impact  parameters  and  the  response  parameters.  This  re- 
lationship can  be  used  to  make  design  tradeoffs  between  the  penetrator 
structure  and  the  Impact  conditions.  For  example,  for  a given  penetrator 
and  e given  target,  a tradeoff  can  be  made  between  Impact  velocity  and 
angle  of  attack. 
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This  Impact-response  relationship,  including  both  terradynamics 
and  structural  response,  is  complex  and  not  understood  for  many  impact 
conditions.  However,  some  insight  can  be  gained  by  introducing  an 
appropriate  set  of  parameters  that  describe  the  loading  on  the  pene- 
trator.  The  impact-load  relationship  (determined  only  from  terradynamics) 
and  the  load-response  relationship  (determined  entirely  by  structural 
response)  can  be  treated  separately  and  combined  to  produce  the  impact- 
response  relationship  sought. 

The  load-damage  relationship  can  be  determined  readily  using  the 
analysis  described  above;  such  a relationship  for  a simple  tube  was 
shown  in  Figure  6.  In  contrast,  the  terradynamic  impact-load  relation- 
ship is  not  as  well  understood.  Therefore,  we  make  two  reasonable 
assumptions,  based  on  currently  available  Information,  to  obtain  an 
impact-load  relationship  that  allows  this  procedure  to  be  Illustrated 
and  shows  the  nature  of  the  impact-response  relationship.  From  available 
theoretical  and  experimental  results,  we  first  make  the  approximation 
that  the  axial  loading  force  F la  proportional  to  the  Impact  velocity 
V.  For  normal  impacts  this  relationship  has  been  verified  for  soil 

targets  and  also  appears  to  hold  for  rock  targets.*^  Second,  we  assume 

+ 

that  the  lateral  load  is  proportional  to  the  angle  of  attack  a, 
over  the  range  of  a of  interest.  This  assumption  is  consistent  with 
a comparison  between  predicted  strain  and  strain  measured  in  angle-of- 
attack  reverse  ballistics  tests  performed  by  AVCO. 

With  this  impact-load  relationship  we  transform  the  load-response 
relationship  into  the  critical  Impact  curves  shown  in  Figure  1 for  four 
penetrator  structures,  These  curves  are  plots  of  combinations  of  im- 
pact velocity  V and  angle  of  attack  a for  which  the  peak  compressive 
stress  in  the  penetrator  is  constant.  Thus,  the  curves  give  the  tradeoff 
between  Impact  velocity  and  angle  of  attack.  The  specific  curves  drawn 

Vi  — — - 

P.  F.  Hadala,  private  communication,  January  1977. 
t -1 

More  precisely,  we  assume  that  tan  n - 3a  where  n “ F /F^,  as 
discussed  in  Section  5. 
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assume  that  the  critical  stress  o Is  ten  times  the  "pressure" 

max 

^1000  " ^1000^^  ® m/sec)  normal  impact.  Structure  B 

Is  similar  to  that  already  used  for  deep  penetration  Into  hard  targets. 
Structure  C is  similar  to  that  proposed  for  shallow  penetration. 
Structures  A and  D were  analyzed  to  show  the  effects  of  changing  f,/d 
or  a/h. 

The  critical  impact  curves  can  'be  used  to  select  the  penetrator 
material,  on  the  basis  of  yield  stress,  for  a given  range  of  V and  a 
within  which  the  penetrator  must  function.  For  example,  if  for 
Structure  B and  a given  target,  the  maximum  impact  velocity  is  to  be 
2000  ft/sec  and  the  maximum  angle  of  attack  is  to  be  3 degrees,  a 
material  with  a yield  strength  of  at  least  10  is  needed.  In  a 

similar  way  the  critical  load  curves  could  be  used  to  select  targets 
for  which  a given  system  (i.e.,  specified  o^,  V,  and  «)  could  be  used, 

Critical  Impact  curves  can  also  be  used  to  interpret  results  of 
experiments  of  more  detailed  load  and  structural  response  calculations 
(e.g.  , a finite  element  code  prediction).  A particular  experiment  or 
detailed  calculation  gives  a single  point  on  a critical  impact  curve. 
Although  such  a data  point  la  probably  more  accurate  than  the  curves 
calculated  with  the  beam-mass  model  used  here,  more  points  are  needed 
to  determine  the  shape  of  the  curve.  The  curves  calculated  with  the 
beam-mass  model, can  be  used  to  determine  this  shape.  Thus  the  appro- 
priate values  of  V and  u can  be  selected  to  minimize  the  data  points 
needed  to  define  the  curve.  For  example,  for  Structure  B small  incre- 
ments In  (X  and  larger  increments  in  V should  be  made  to  efficiently 

define  the  a /Pm^n  " 10  curve  for  small  a. 
max  1000 

The  critical  load  curves  are  perhaps  most  useful  for  making  design 
tradeoffs  among  candidate  penetrator  structures.  For  example,  the 
curves  for  Structure  B (deep  penetrator)  and  for  Structure  C (Shallow 
penetrator)  pass  through  a similar  region  in  the  V-u  plane  and  Inter- 
sect at  V ■ 1500  ft/sec  and  u ■ 5.7  degrees.  However,  the  curve  for 


Structure  B is  steeper  and,  for  small  angles  of  attack,  this  structure 
can  withstand  greater  Impact  velocities.  Both  these  penetrators  can  be 
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made  stronger  by  Increasing  only  the  wall  thickness  (from  that  of 
Structure  C to  that  of  Structure  A)  or  decreasing  the  length  (from  that 
of  Structure  B to  that  of  Structure  A).  Note  that  either  of  these 
changes  reduces  the  payload  volume.  Also,  for  small  angles  of  attack, 
the  velocity  range  can  be  extended  more  by  Increasing  wall  thickness 
than  by  decreasing  length,  Similar  conclusions  can  be  drawn  by  comparing 
curves  for  Structures  B and  C with  that  of  Structure  D. 

The  method  developed  here  for  characterizing  the  strength  of 
penetrator  structures  has  been  applied  to  four  idealized  structures. 

The  principal  advantage  of  this  procedure  is  that  it  characterizes  the 
response  over  a range  of  Impact  conditions.  It  can  be  applied  to  more 
complex  structures  than  those  analyzed  here  to  plan  and  Interpret 
experiments  and  detailed  calculations  and  to  compare  the  performance  of 
different  structures  for  a particular  application, 

REPORT  ORGANIZATION 

The  remainder  of  this  report  is  divided  into  four  sections. 

Section  3 describes  the  experiments  on  1/4-soale  model  penetrator-llke 
structures  end  the  measured  response.  Section  4 presents  the  develop- 
ment of  the  mathematical  analysis,  comparison  of  predicted  response 
with  measured  response,  and  the  effect  of  certain  parameters,  such  as 
loading  rise  time,  on  the  predicted  response.  Section  5 presents 
critical  Impact  curves,  that  is,  combinations  of  Impact  velocity  and 
angle  of  attack  that  produce  the  same  peak  stress  in  a given  structure. 
Structures  analyzed  have  radlus-to-thickness  ratios  ranging  from  2.0  to 
10.0  and  length-to-diameter  ratios  ranging  from  3.6  to  8.0. 
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3.  SIMULATOR  EXPERIMENTS 

The  primary  objective  of  the  experimental  work  described  here  was 
to  evaluate  a concept  for  economical  testing  of  penetrator  structures. 

To  make  such  an  evaluation  the  test  fixture  was  designed  for  loading 
nominal  1/4-ncale  model  penetratora;  the  resulting  design  can  be 
fabricated  In  full  scale  for  testing  actual  penetrator  structures.  A 
secondary  objective  was  to  provide  experimental  data  for  comparison  with 
the  analysis. 

The  load  simulator  design  Is  based  on  two  basic  principles.  First, 
damage  to  the  penetrator  occurs  during  or  shortly  after  the  load  reaches 
its  peak  amplitude;  therefore,  only  this  portion  of  the  load  history 
need  be  simulated.  The  change  In  momentum  of  the  penetrator  during  this 
time  is  small  compared  to  the  total  momentum  of  a penetrator  on  impact. 
Therefore  by  applying  only  this  early  time  portion  of  the  load  to  a 
penetrator  initially  at  rest,  the  response  of  Interest  is  simulated  but 
the  resulting  momentum  transmitted  to  the  structure  Is  much  less  than 
the  initial  momentum  required  in  ballistic  or  reverse  ballistic  tests. 

This  allows  simulator  tests  to  be  performed  in  the  laboratory.  Second, 
observed  maximum  compressive  strains  and  failures  In  penetrator  structures 
occur  at  distances  greater  than  one  diameter  from  the  penetrator  nose. 

The  stress  distribution  over  the  cross  section  at  these  locations 
depends  only  on  the  resultant  forces  applied  to  the  end  of  the  structure. 
Thus  the  load  simulator  was  designed  to  produce  the  resultant  axial  and 
lateral  loading  forces  with  the  appropriate  time-history,  and  not  the 
details  of  the  pressure  load  distribution  on  the  penetrator  nose. 

The  remainder  of  this  section  covers  three  aspects  of  the  load 
simulation  tests;  the  design  and  construction  of  the  test  fixture,  the 
1/4-scale  model  structures  that  were  tested,  and  some  typical  results  of 
the  teste. 


I 
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TEST  FIXTURE 


Operation 

Figure  7 Is  a sectioned  assembly  drawing  of  the  fixture  In  con- 
figuration for  simulating  normal  Impact  loads  on  1/4-scale  model  pene- 
trators . The  fixture  operates  as  follows:  high-pressure  gaseous 
explosive  products  are  produced  in  the  explosive  chamber  by  detonation 
of  a solid  explosive.  The  gas  flows  through  an  orifice  plate  and  Into 
a cylinder  containing  a piston  that  Is  In  contact  with  the  penetrator. 

The  load  Is  transmitted  to  the  penetrator  through  the  piston  (the  details 
of  this  Interface  are  discussed  later).  The  fixture  was  designed  so  that 
the  rise  time  of  the  pressure  load  could  be  varied  by  using  different 
Initial  piston  displacements  or  different  orifice  areas.  The  duration 
of  the  nearly  constant  load  plateau  can  be  varied  by  using  vent  holes  at 
different  locations  along  the  cylinder,  and  the  decay  time  can  be  varied 
by  using  different  size  vent  holes. 

To  aid  In  the  design  of  the  fixture,  the  pulse  produced  by  a given 
geometry  was  predicted  by  using  the  GASLEAK  computer  code,®  which  models 
the  flow  of  gases  in  a series  of  chambers  connected  by  orifices.  The 
theoretical  model  of  the  flow  assumes  that  (1)  the  duration  of  the 
loading  pulse  is  long  compared  with  the  transit  time  of  pressure  waves 
In  each  chamber  (quasi-steady  flow),  (2)  neglit;iblc  heat  is  transferred 
from  the  hot  gas  to  the  surrounding  cylinder  (adiabatic  flow) , and 
(3)  the  hot  detonation  products  behavo  as  a perfect  gas.  Experience  in 
other  similar  applications  has  shown  excellent  agreement  between  the 
theor telcally  predicted  pressure  pulse  and  the  experimentally  measured 
pulse. 

Construction 

The  device  is  constructed  on  a stack  of  alloy  steel  rings  and 
circular  plates  clamped  together  by  eight  tie  rods.  The  explosive 
chamber  Is  formed  by  a thick  annulus  that  fits  between  the  base  plate 
and  the  orifice  plate.  The  cross-sectional  area  of  each  of  the  six 
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orifices  may  be  varied  by  Inserting  plugs  drilled  with  the  proper  size 
hole  Into  the  recesses  on  the  lower  side  of  the  orifice  plate.  The 
cylinder  In  which  the  piston  travels  is  formed  by  the  vent  ring  and 
the  cylinder  plate.  The  Initial  volume  of  the  piston  chamber  is 
controlled  by  a spacing  washer  between  the  piston  and  orifice  plate, 
which  sets  the  initial  piston  displacement. 

The  vent  ring  has  three  pairs  of  vent  holes,  each  pair  at  a 
different  axial  location.  These  holes  vent  the  piston  chamber  to  the 
atmosphere  as  the  piston  passes.  The  size  of  the  vent  area  may  be 
varied  with  threaded  plugs  that  reduce  the  area  of  the  vent  holes  or 
close  the  holes  completely.  The  axial  location  of  the  holes  Is  set 
by  the  thickness  of  the  spacing  ring.  Two  pressure  gages  are  mounted 
diametrically  opposed  In  the  vent  ring  to  measure  the  chamber  pressure. 
After  the  pulse  Is  produced,  the  piston  decelerates  by  Impacting  the 
energy-absorbing  aluminum  honeycomb. 

Experimental  Setup 

To  simulate  the  load  on  a penetrator  that  Impacts  at  an  angle  of 

attack,  the  piston  Is  designed  to  produce  s combined  axial  and  lateral 

loading  and  to  measure  the  load  applied  to  the  penetrator.  The  combined 

loading  Is  produced  by  tilting  the  penetrator  through  an  angle  6 with 

respect  to  the  piston  axis,  as  shown  in  the  sectioned  drawing  of  the 

piston  in  Figure  8.  The  angle  0 may  be  varied  from  zero  to  20® 

(0.35  rad).  As  discussed  In  Section  3,  this  corresponds  to  an  angle 

of  attack  of  about  3.7  degrees.  In  structural  response  tests,  the 

* 

vertical  force  Is  measured  with  a piezoelectric  load  cell 

(Kistler  906A) . Two  discs  arc  placed  between  the  penetrator  and  the 
load  cell,  and  the  Interface  between  the  discs  is  lubricated  with  a 


As  discussed  later,  a load  cell  was  not  used  in  the  long-rod  calibration 
tests , 
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FIGURE  8 PISTON  FOR  ANGULAR  LOADING 
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high-pressure  solid  lubricant  (Dow  Corning  Molykote  321R)  to  minimize 
the  horizontal  frictional  force.  When  the  load  Is  applied,  sliding 
occurs  between  the.  two  discs,  and  the  vertical  force  and  a hori- 
zontal frictional  force  are  applied  to  the  penetrator, 

The  ratio  of  the  applied  loads  be  determined  from  the 

motion  of  the  penetrator.  With  the  assumption  that  this  ratio  Is 

+ 

constant  during  application  of  the  loading,  the  load  ratio  Is  related 
to  the  angle  0 between  the  path  of  the  center  of  mass  of  the  pene- 
trator and  the  vertical  axis  of  the  simulator  as  follows: 

F 

- tan  3 

The  angle  3 was  measured  using  high  speed  photography;  the  largest 
value  of  3 measured  was  about  9 degrees.  Figure  9 shows  the  assembled 
device  with  a 1/4-scale  model  penetrator  In  position  for  a normal  Impact 
simulation.  Before  the  experiment,  the  model  Is  held  In  place  by  low- 
strength,  machined  Styrofoam  rings  that  fit  Inside  the  top  end  plate 
and  the  cylinder  plate.  In  angular  Impact  tests,  the  Initial  angle  6 
between  the  axis  of  the  penetrator  and  the  axis  of  the  simulator  is 
measured  with  a vernier  protractor.  After  the  load  simulation,  the 
penetrator  leaves  the  device  and  Is  stopped  by  an  external  energy 
absorber  (aluminum  honeycomb  or  Styrofoam)  located  In  a 5-foot-long 
(1. 52-m-long)  safety  shroud.  The  shroud  ensures  containment  of  the 
model  penetrator  after  the  simulation.  Figure  10  shows  the  device  with 
the  safety  shroud  In  place  for  testing. 


For  Impact  at  an  angle  of  attack,  the  component  of  the  friction  force 
lateral  to  the  penetrator  acts  in  the  opposite  direction  from  the 
lateral  force  we  wish  to  simulate. 


i* 

This  assumption  was  examined  In  a series  of  tests  In  which  the  load 
duration  was  shortened  by  means  of  a mechanical  stop  placed  above  the 
piston.  The  tests  Indicate  that  Che  load  ratio  variation  during  the 
Impact  simulation  is  within  measurement  error. 
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FIGURE  9 ASSEMBLED  TEST  FIXTURE  FOR  SIMULATING  IMPACT  LOADS 


In  some  experiments  high  speed  movies  were  taken  to  measure  the 
angle  0.  In  these  experiments  the  camera  (Hycam  Model  410004)  was 
placed  3 feet  (91  cm)  from  the  simulator  and  photographed  12  inches 
(30  cm)  of  penetrator  travel  at  the  nominal  rate  of  lOiOOO  frames  per 
second  (a  safety  shroud  that  allows  a large  field  of  view  was  used 
rather  than  the  shroud  shown  in  Figure  10).  The  center  of  mass  of  the 
penetrator  was  marked  with  a bullseye  to  allow  its  path  to  be  photo- 
graphed after  the  load  simulation.  To  provide  a reference  for  tracing 
the  motion  of  the  model  penetrator,  a grid  was  placed  on  the  top  of 
the  simulator,  filmed,  and  removed  before  the  experiment. 

After  the  experiment  the  motion  of  the  center  of  mass  of  the  model 
was  determined  by  using  a Telereadex  film  analyzer.  A frame  of  the  film 
showing  the  reference  grid  was  first  projected  onto  a table  and  traced 
onto  a sheet  of  paper.  Subsequent  frames  were  projected  onto  the  traced 
grid,  The  position  of  the  center  of  mass  in  several  frames  was  marked 
on  the  grid.  Typically,  the  path  of  motion  was  about  15  inches  (38  cm) 
long  on  the  projection  screen.  The  tangent  of  the  angle  0 between 
the  path  of  the  center  of  mass  and  a vertical  grid  line  was  then 
measured  directly.  Accuracy  of  measurement  was  about  ±0,01  radians, 

MODEL  STRUCTURES 

Two  general  types  of  structures  were  used  in  the  experiments. 

First,  we  used  a calibration  rod  long  enough  that  waves  reflected  from 
its  free  end  did  not  reach  the  piston  during  the  load  rise  time.  Then, 
we  tested  several  shorter  structures  typical  of  penetrators. 

Calibration  Rod 

In  the  calibration  experlmenCs  (Tests  12  through  30)  we  used  an 
AISI-1020  steel  rod  1.5  inches  (3.81  cm)  in  diameter  and  30  inches 
(76.2  cm)  long.  To  simulate  a normal  impact,  we  placed  this  long  rod 
perpendicular  to  the  top  face  of  the  piston  with  the  rod-piston  interface 
conditions  shown  in  Figure  11(a).  The  3/4-'inch-diameter  (1. 90-cm-diameter) 
steel  disc  applied  the  load  to  the  central  area  of  the  rod  so  that  any 
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FIGURE  11  LOADING  CONFIGURATION  FOR  LONG  ROD 
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eccentricity  of  the  loading  would  be  reduced.  The  lead  sheet  smoothed 
any  surface  irregularities.  Strain  measured  near  the  front  end  of  the 
long  rod  was  used  as  a check  on  the  load. 

The  calibration  rod  Is  Instrumented  with  eight  strain  gages 
(Micro-Measurements  Type  EP-08-250GB-120)  to  measure  the  rod  strain  In 
the  axial  direction.  The  locations  of  the  gages  are  shown  In  Figure  11(b). 
The  four  gages  at  each  axial  station  are  evenly  spaced  around  the  rod. 


Penetrator  Models 

Figure  12  Is  a general  schematic  drawing  of  four  of  the  penetrator 
models  used  In  the  experiments.  The  models  have  a 1.50-lnch  (3.81-cm) 
outer  diameter  and  have  a hemispherical  front  end  rather  than  a pointed 
nose  to  allow  the  desired  load  to  be  applied  more  easily  and  accurately. 
Each  model  Is  instrumented  with  eight  strain  gages  to  measure  strain  in 
the  axial  direction.  The  details  of  each  model  are  given  In  Table  1. 

The  solid  steel  model  is  9 Inches  (22.9  cm)  long  and  Is  made  of  AISI 
1020  steel.  It  Is  the  simplest  model  since  it  has  a uniform  cross  section 
(except  at  the  loaded  hemispherical  end) . The  solid  aluminum  model  Is 
geometrically  Identical  to  the  solid  steel  model  but  is  made  of  6061-T6 
aluminum.  It  was  only  used  to  check  the  behavior  of  the  piston  sliding 
plate  when  a lightweight  model  was  tested;. 

i 

Two  of  the  models  have  cavities  (indicated  by  the  dashed  lines  in 
Figure  12).  Both  are  made  of  AISI  1020  steel.  Ihe  thick-walled  model 
Is  9 Inches  (22.9  cm)  long  with  a 0 .875-lnch-diameter  (2. 22-cm-dlameter) 
and  a 6-inch-long  (15.24-cm-long)  cylindrical  cavity.  The  radlus-to- 
thickness  ratio  of  the  cylindrical  portion  of  this  model  is  typical  of 
deep  earth  penetrator  structures.  The  thin-walled  model  is  6.625  Inches 
(16.8  cm)  long  with  a 1. 305-inch-diemeter  (3. 43-cm-diameter)  and  a 
4. 7-lnch“long  (11. 94-cm-long)  cavity.  The  radlus-to-thickness  ratio  of 
this  model  Is  typical  nf  proposed  shallow  penetrator  structures. 
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On  the  solid  steel  model  the  positions  of  gages  2 and  6 are  Interchanged 
with  those  of  gages  4 and  8. 

MA-3091-iaiA 


FIGURE  12  MODEL  PENETRATOR  WITH  CAVITY 


Table  1 

MODEL  PENETRATOR  STRUCTURES 


Strain  Gage 
Cavity  Louatlone^ 


Model 

Abbr. 

Mass 

(gm) 

Length 

Length 

(cm) 

Diam . 
[I.D.  ] 
(cm) 

Front 

[Sp] 

(cm) 

Aft 

tSA] 

issil 

Solid  1020  steel 

SS 

2030 

22.86 

— 

“ 

3.81 

11.43 

Thick-walled 

1020  steel 

TkW 

1535 

22.86 

15.24 

2.22 

3.81 

11.43 

Thin-walled 

1020  steel 

TnW 

620 

16.83 

11.94 

3.43 

6.10 

10.45 

Solid  6061-T6 
aluminum 

SAl 

683 

22.86 

__ 

— 

— 

Tapered  6061-T6 

TAl 

759 

19.05 

__ 

5.08 

9.53 

aluminum 


Letters  In  brackets  refer  to  achewatlc  drawing  of  model  penetrator  In 
Figure  12. 

^Distance  measured  from  front  of  model  structure. 
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In  four  tests  a tapered  solid  aluminum  model  was  used  to  provide 
data  for  future  analysis  of  tapered  structures.  The  loaded  end  Is  a 
1.50-inch-dlameter  (3. 81~cm-dlameter)  hemisphere  similar  to  the  other 
models.  The  body  tapers  outward  to  1.875-lnch-diameter  (4. 76-ctn-diameter) 
at  the  aft  end.  The  overall  length  is  7.50  inches  (19.0  cm). 

EXPERIMENTAL  RESULTS 

The  details  of  each  test  and  the  responses  measured  are  given  In 
Appendix  A (Tables  A-1  through  A-4) . The  strain  responses  measured  In 
Tests  64  and  66  are  compared  with  the  predicted  responses  in  Section  4. 

The  test  results  reflect  the  development  and  refinement  of  the 
simulation  technique.  For  examplet  Initial  tests  with  the  long  rod 
indicated  that  any  slight  deviation  in  the  perpendicularity  of  the  rod 
and  the  piston  face  caused  the  load  to  be  applied  eccentrically  at  the 
edge  of  the  rod.  This  effect  varied  from  test  to  test,  causing  irregular 
bending  strains.  A small-diameter  interface  disc  was  placed  bem^een  the 
rod  and  the  piston  to  reduce  this  effect.  Although  this  reduced  the 
bending,  significant  undesired  bending  strains  still  occurred.  The 
problem  was  resolved  in  the  later  models  by  using  a hemispherlcfal  loaded 
end.  Other  parts  of  the  simulation  technique  added  during  the  test  pro- 
gram included:  a load  cell  for  measuring  the  applied  vertical  load  and 
high  speed  photography  for  determining  the  ratio  of  the  applied  hori- 
zontal force  to  the  applied  vertical  force. 

In  the  tests  with  angular  loading,  two  opposing  strain  gages  at 
each  of  the  two  stations  were  aligned  in  the  plane  of  bending.  The 
other  two  gages  at  each  station  were  in  the  neutral  plane  in  which  no 
bending  occurs.  In  the  figures  that  follow,  the  four  gages  from  each 
station  are  grouped  together  as  indicated  in  Figure  4 : the  strain 
gage  record  from  each  station  that  contains  the  compressive  side  of 
bending  is  displayed  in  the  upper  right-hand  corner  of  the  group  of 
four  records  from  that  station.  The  gage  record  containing  the  tensile 
side  of  bending  is  in  the  lower  left-hand  corner  of  the  group.  The 
other  two  records  in  each  group  show  only  the  axial  strain.  (Note 
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that  the  gage  numbers  of  the  gages  In  the  bending  plane  varied  from 
test  to  test.)  When  all  eight  strain  gage  records  are  shown,  the  front 
station  Is  the  top  group  of  four,  and  the  aft  station  Is  the  lower  group. 

Long  Rod 

Figure  13  shows  the  chamber  pressure  and  front  station  strain 
records  from  Test  16.  In  this  test  the  long  calibration  rod  was  loaded 
axially.  The  nonuniformity  of  the  strains  due  to  the  eccentrically 
applied  load  discussed  previously  can  be  seen  by  comparing  the  magnitudes 
and  shapes  of  the  four  strain  records.  The  average  peak  strain  of  0.020 
percent  corresponds  to  an  axial  force  of  10,600  pounds  (47.1  kN).  The 
peak  pressure  corresponds  to  an  axial  force  of  11,300  pounds  (50.3  kN) . 

The  difference  In  the  two  force  levels  is  caused  by  the  inertia  of  the 
piston. 

Solid  Steel  Model 

Figure  14  shows  the  pressure,  load  cell,  and  strain  gage  data 
recorded  In  Test  37  In  which  the  9-lnch-long  (22 .9-cm-long)  solid  steel 
model  was  normal  to  the  piston  face.  The  results  of  this  test  show 
several  characteristics  typical  of  the  model  tests.  As  the  load  is 
applied  to  the  model,  the  hemispherical  end  Indents  the  sliding  disc 
on  the  piston,  causing  plastic  deformation  of  the  disc.  This  deformation 
and  the  inertia  of  the  piston  account  for  the  difference  betv7een  the 
rise  time  shown  on  the  pressure  record  and  that  shown  on  the  load  cell 
record.  The  load  cell  record  shows  that  the  load  rises,  gradually  levels, 
and  then  drops  off.  The  simulation  load  comprises  the  load  rise  and 
leveling  off  and  ends  as  the  load  begins  to  drop.  As  will  be  discussed 
in  Section  4,  the  peak  penetrator  response  occurs  shortly  after  the  load 
reaches  its  peak;  therefore,  the  duration  of  the  simulation  is  adequate. 
The  uniformity  of  the  axial  strain  at  each  station  is  attributed  to  the 
hemispherical  end.  This  test  and  similar  ones  show  that  the  circum- 
ferential variation  in  measured  strain  due  to  axial  loads  is  about 
3 percent. 
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400  psi/cm  (2758  kPa/cm)  200  /usuc/cm  400  psi/cm  (2758  kPa/cm)  50  /usec/cm 

(a)  PRESSURE  RECORD  (GAGE  1)  (b)  PRESSURE  RECORD  (GAGE  1) 


GAGE  2 


GAGE  3 


0,0079%/cm  50  psec/cm  0.0079%/cm  60  /jsec/cin 

(c)  STRAIN  GAGE  RECORDS  FROM  FRONT  STATION 

MP-3091-118AI 


FIGURE  13  TEST  16— CALIBRATION  TEST  WITH  LONG  ROD 
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1000  psi/cm  (6894  k Pa/cm) 

60  /isec/cm  100  psec  delay 
(a)  PRESSURE  RECORD  (GAGE  2) 


5000  Ibs/cm  (22.2  kN/cm) 
50  psec/cm  200  /usee  delay 

(b)  LOAD  CELL  RECORD 


-H, 


tiAUfl 

■ 

(c)  0.01  %/cm  60  /usec/cm  200  /usee  DELAY 

STRAIN  GAGE  RECORDS  FROM  FRONT  STATION 


I 0.01  %/cm  60  /usec/cm  200  psec  delay 

(d)  STRAIN  GAGE  RECORDS  FROM  AFT  STATION 

1 IVtA-3091-136A 

I FIGURE  14  TEST  37— NORMAL  LOADING  OF  SOLID  STEEL  MODEL 
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Figure  15  shows  the  results  of  Test  34.  In  this  test  the  axis  of 
the  solid  steel  model  was  tilted  5 degrees  (0.09  rad)  with  respect  to 
the  axis  of  the  piston,  and  the  explosive  charge  mass  was  the  same  as 
that  of  Test  37.  The  bending  strain  due  Co  the  resulting  transverse 
load  can  be  seen  most  clearly  In  the  records  from  Gages  6 and  8 at  the 
aft  station.  Gages  5 and  7 were  on  the  neutral  plane  of  bending  and 
therefore  recorded  only  the  axinl  strain  of  the  rod{  the  peak  axial 
strain  of  0.015  percent  is  comparable  to  Chat  for  normal  loading 
(Figure  15) . The  record  from  Gage  8 shows  chat  In  this  case  the  super- 
posed compressive  axial  strain  and  the  tensile  bending  strain  almost 
cancelled  each  other,  resulting  In  a small  total  strain  during  the  500 
Visec  recorded.  Thus  the  peak  bending  strain  was  also  about  0.015  percent, 

Thick-Walled  Model 

Figure  16  shows  the  applied  load  and  aft  station  strain  records 
from  Test  40  in  which  the  3.37-pound  (1.53-kg)  thiek-wa.Ued  steel  model 
was  normal  to  the  piston.  The  peak  load  was  16,500  pounds  (73.4kH)  and 
Che  average  axial  strain  was  0.027  percent.  As  in  Test  37,  the  cir- 
cumferential variation  In  strain  Is  small. 

Figure  17  and  18  show  the  results  from  Tests  43  and  64.  In  Test  43 
the  peak  load  was  16,000  pounds  (71.2  kN)  and  the  peek  axial  strain  at 
the  aft  station  was  0.023  percent.  In  Test  64  the  peak  load  was 
15,000  pounds  (66.7  kN)  and  the  peak  axial  strain  was  0,022  percent. 

The  initial  angle  of  tilt  6 In  Test  64  tlS^SO'  (0.340  rad)]  was  2.7 
times  that  of  Test  43  [7® 10*  (0.125  rad)].  The  peak  bending  strain  in 
Test  64  (0.042  percent)  was  2.6  times  that  of  Test  43  (0.016  percent). 

Thus  there  Is  a nearly  linear  relation  betwaen  the  angle  of  tilt  C>  and 
the  bending  strain. 
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1000  psi/cm  (6894  kPa/cm) 

50  A<sec/cm  100  fxsec  delay 
(a)  PRESSURE  RECORD  (GAGE  1) 


X, 

/ 1 

Fraw 

''  Mf': 


6000  Ibs/cm  (22.2  kN/cm) 
50  Msec/cm  150  /usec  delay 
(b)  LOAD  CELL  RECORD 


(c)  0.01  %/cm  50  /Liseo/cm  160  /asac  DELAY 

STRAIN  GAGE  RECORDS  FROM  FRONT  STATION 


J I 


CIBMJ 


■ — . 

IHMiMI 

\ 

FIGURE  15 


Enijaij 


0,01%/cm  50  (jsec/cm  160  /usee  delay 
(d)  STRAIN  GAGE  RECORDS  FROM  AFT  STATION 

MA-3091  -124A 

TEST  34-5  (0.09  radl  ANGULAR  LOADING  OF  SOLID  STEEL 
MODEL 
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.LOAD  CELL  RECORD 
[5000  Ibs/cm  (22,2  kN/cm) 
50 /jsec/cm  150 /j sec  delay] 


GAGE  8 


GAGE  5 


(b)  STRAIN  GAGE  RECORDS  FROM  AFT  STATION 
(0.016%/cm,  50  /usec/cm,  150  ;usec  delay) 

MA-3Qn1-i;K]A 

FIGURE  16  TEST  40-  NORMAL  LOADING  OF  THICK-WALLED  MODEL 


(a)  LOAD  CELL  RECORD 

[5000  Ibs/cm  (22,2  kN/cin) 
50  iDSec/cm  150  jusec  delay] 


GAGE  7 GAGE  6 


GAGE  8 GAGE  5 

(b)  STRAIN  GAGE  RECORDS  FROM  AFT  STATION 
(0,016%/cm,  50  /Jsec/cm,  150  ;jscc  cleltiy) 


MP-;)0Oi-ia;'A 

FIGURE  17  TEST  43-7''10'  (0.125  rad)  ANGULAR  LOADING  OF 
THICK-WALLED  MODEL 
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Thin-Walled  Model 


Figures  19  and  20  show  the  records  from  Tests  55  and  66  in  which 
the  1.37-pound  (0.62-kg)  thin-walled  model  was  loaded  both  axially 
and  laterally.  In  teat  55  the  angle  of  tilt  was  8®10'  (0.143  rad), 
the  peak  load  was  11,000  pounds  (48.9  kN),  the  peak  axial  strain  was 
0.019  percent,  and  the  bending  strain  was  small  relative  to  the  axial 
strain.  In  Test  66  the  angle  of  tilt  was  18“15'  (0.319  rad),  the  peak 
load  was  11,250  pounds  (50.0  kN) , the  peak  axial  strain  was  0,023  percent, 
and  the  peak  b.;ndlng  strain  was  0.015  percent.  (Note  that  the  strain 
records  have  different  time  scales  in  the  two  tests.) 

CONCLUSIONS 

The  test  results  show  that  the  load  simulator  can  apply  to  model 
penetratora  loids  similar  to  those  occurring  in  angle-of-attack  impacts, 
The  structural  response  of  the  penetrator  has  been  accurately  measured 
using  hard-wired  strain  gages.  The  tests  indicate  that  the  loader 
could  also  be  built  in  a larger  size  to  test  full-scale  penetrators. 

For  the  1.5-inch-diameter  penetrator  models  tested  here,  the 
shortest  rise  time  produced  was  about  150  ysec.  For  a penotrator  whose 
nose  length  is  about  1,5  times  its  diameter,  this  corresponds  to  an 
impact  velocity  of  1.5  (1.5  inch)/l50  Usee  » 1250  ft/sec  (385  m/sec). 

It  may  be  desirable  to  further  decrease  the.  load  rise  time  by  making 
the  orifice  plate  thinner  or  by  making  the  piston  surface  harder  to 
simulate  greater  Impact  vclocoties. 
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(a)  LOAD  CELL  RECORD 

[2500  Ibs/cm  (11.1  kN/cm) 
50  psec/cm  150  ^sec  delay] 


GAGE  6 


GAGE  5 


GAGE  7 GAGE  8 

(b)  STRAIN  GAGE  RECORDS  FROM  AFT  STATION 
(0.016%/cm,  50  psec/cm,  150  i^soc  delay) 

MP-30ni-139 


FIGURE  20  TEST  ee-IB-^IB’  (0,319  rad)  ANGULAR  LOADING  OF 
THIN-WALLED  MODEL 
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4 ANALYSIS  OF  PENETRATOR  RESPONSE 


In  this  section  we  develop  the  theory  for  predicting  penetrator 
response  under  angle-o£-attack  Impacts » compare  the  response  predicted 
by  this  theory  with  that  measured  In  simulation  experiments  and  with 
that  measured  in  a reverse  ballistics  test  performed  at  AVCO,  and 
investigate  the  effect  of  specific  loading  and  structural  parameters 
on  penetrator  response, 

As  discussed  previously,  large  strains  and  failures  in  penetrators 
were  observed  experimentally  In  several  angle-of-attack  impacts. 
Comparison  of  the  location  of  these  failures  with  the  peak  strains 
measured  in  the  simulator  tests  performed  here  and  in  reverse  ballistics 
tests  performed  at  AVCO  indicate  that  these  failures  occur  as  follows. 
Impacts  at  an  angle  of  attack  produce  bending  as  well  as  axial  thrust 
in  the  penetrator.  When  the  total  compressive  stress  due  to  axial  thrust 
plus  bending  reaches  the  yield  stress,  the  bending  stiffness  of  the 
penetrator  is  reduced  and  the  bending  strain  increases  further  until 
failure  occurs.  Thus,  the  analysis  is  based  on  elastic  response  because 
that  is  what  determines  whether  or  not  fallu-re  can  ultimately  occur. 

THEORY 

The  central  portion  of  the  penetrator,  drawn  schematically  in 
Figure  21(a),  is  modeled  by  elastic  Timoshenko  beam  theory  including 
thrust.  This  formulation  includes  the  two  dominant  response  mechanisms 
governing  angle-of-attack  impacts!  axial  compression  and  bending  under 
lateral  loads.  The  front  portion  of  the  solid  nose  and  the  solid  aft 
end  are  modeled  as  rigid  masses.  Loading  of  this  structure  is  specified 
by  the  magnitude,  time  history,  and  orientation  of  a resultant  impact 
force;  in  all  cases  studied  here  we  take  the  resultant  force  to  pass 
through  a fixed  point  in  the  front  mass. 


54 


(b) 

FIGURE 


FORCES  AND  MOMENTS  ON  THE  BEAM  AND  FRONT  AND  AFT 
RIGID  MASSES 

MA-309i-109A3 


21  MATHEMATICAL  MODEL  OF  AN  IDEALIZED  PENETRATOR 
STRUCTURE 
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Beam  Equations 


Six  equations  governing  the  motion  of  the  beam  section  are  re- 
quired; two  governing  axial  compression  and  four  governing  bending. 
The  equations  governing  axial  response  of  the  beam  are 

N - EA  - 0 (la) 


3x 


P A 


3t2 


(lb) 


where  N Is  axial  force  In  compression,  E Is  modulus  of  elasticity, 
A Is  cross-sectional  area,  u Is  axial  displacement,  x is  the 
coordinate  along  the  beam  axis,  p Is  density,  and  t Is  time.  The 
first  of  these  equations  Is  Hooke's  law  and  the  second  Is  the  equation 
of  motion  In  the  axial  direction. 


The  equations  governing  bending  response  of  the  beam  are 


M + El  - 

0 

(Ic) 

/AG(yj^  - f)  - 

0 

(Id) 

» + ’ J ■ 

0 

(le) 

ia  , DA  - 

0 

(If) 

where  M is  bending  moment,  I is  cross-sectional  moment  of  inertia, 
V is  that  part  of  the  slope  of  the  deflection  curve  caused  by  bending 
(the  total  slope  y^  is  caused  by  both  bending  and  shearing),  Q Is 
shear  force,  k'  - ir^/12  Is  the  shear  correction  coefficient,  G is 
shear  modulus,  and  y is  transverse  deflection.  The  first  of  these 
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equations  la  the  bending  moment-curvature  relation  of  elementary  beam 
theory.  The  second  expresses  the  shear-force  shear-deformation  relation, 
The  third  is  the  equation  of  motion  for  rotation,  and  the  fourth  is  the 
equation  of  motion  for  translation. 

The  motion  of  the  forward  and  aft  rigid  masses  is  coupled  with 
the  motion  of  the  beam.  It  is  convenient  to  treat  this  interaction  by 
the  characteristic  equations  of  the  beam.  Therefore,  we  now  discuss 
the  characteristic  form  of  equations  (1), 

Equations  (1)  are  a set  of  partial  differential  equations  in  x and 
t.  The  method  of  characteristics  affords  a transformation  of  these  equa- 
tions into  total  differential  equations  that  apply  along  particular,  or 
characteristic,  lines  in  the  x-t  plane.  The  details  of  this  trans- 
formation are  given  in  Reference  7.  The  reciprocal  slopes  of  the 
characteristic  lines  correspond  to  wave  speeds  at  which  disturbances 
can  propagate.  For  the  system  of  equations  (1),  disturbances  can  propa- 
gate at  four  possible  velocities  along  the  beam:  ±c.  and  ±r  , where 

^ D 8 

c.  ■ /b7?  is  the  bur  velocity,  and  c.  ■ \/k''G/p*  is  the  shear  velocity. 

D B 

Axial  waves  propagate  at  ic,  and  are  governed  by  equations  (la)  and 

a 

(lb),  while  transverse  waves  propagate  at  both  ±c,  and  ±c  and  are 

D 8 

governed  by  equations  (Ic)  through  (If)  . To  make  the  subfiequent  analysis 

+ 

more  clear,  we  give  the  label  I to  the  characteristic  lines  having 
slopes  when  referring  to  equations  (la)  and  (lb)  governing  axial 

response,  and  II  to  the  same  lines  when  referring  to  equation  (Ic) 
through  (If)  governing  transverse  response,  us  indicated  in  Figure  22(b). 

4- 

We  give  the  label  III“  to  the  characteristic  lines  having  slopes  ±l/c.  , 

s 

as  also  indicated  in  Figure  22(b).  The  governing  or  churacteristlc  axial 
equations  along  lines  I*  are 


I : d [ N T I ' c A d V ] ■ 0 


where  v ■ 3u/Dt  is  the  axial  particle  velocity. 


Iv., 


+■  + 

The  characteristic  transverse  equations  along  lines  11  ' and  III" 
are 


II*:  d[M  + pi  c^  w]  - ±Cjj  Q dt 
III*:  d[Q  T pA  Cg  W]  - -k'AGw  dt 


(3) 


where  ui  ■ 3'l'/9t  is  the  rotational  velocity  and  w - 3y/dt  is  the 
transverse  velocity. 

Equations  (2)  and  (3)  are  two  independent  sets  of  total  differential 
equations:  that  isi  the  solutions  to  equations  (2)  can  be  found  indepen- 
dently from  the  solutions  to  equation  (3) . Each  set  of  equations  (2) 
or  (3)  are  coupled,  of  course,  and  must  be  solved  as  a group.  The  solution 
to  each  set  was  obtained  by  a forward  differencing  finite  difference  method 
along  the  characteristics.  So  that  solutions  could  be  placed  on  a common 
right  Cartesian  grid,  values  at  points  H and  J in  Figure  22(b)  were 
approximated  from  those  at  points  G and  K by  linear  interpolation 
before  the  numerical  integration,  Thus  the  solutions  at  point  P con 
bo  written  in  terms  of  those  at  points  G and  K.  This  solution  applies 
at  all  interior  points  of  the  grid  of  Figure  22(a).  Details  of  the  solution 
are  given  in  Reference  7.  At  end  points  A and  Z in  the  grid,  the  rigid 
masses  must  be  taken  into  account. 


End  Mass  Equations 

The  coupled  motion  of  the  front  rigid  mass  and  the  front  end  of  the 
Timoshenko  beam  is  governed  by  the  equations  of  motion  and  stress-strain 
relations  of  the  Timoshenko  beam,  by  the  equations  of  motion  of  the  front 
rigid  mass,  and  by  requiring  compatibility  of  displacement  and  rotation  of 
the  front  end  of  the  beam  and  the  front  rigid  mass.  The  equations  for  the 
Timoshenko  beam  are  most  conveniently  represented  by  the  finite  difference 
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(a  COMPLETE  GRID  (COARSE  SPACING  FOR  CLARITY 
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(b)  DETAILS  OF  SHADED  PORTIONS  OF  PART  (a) 

MA-3oei-iBe 

FIGURE  22  SOLUTION  GRID  FOR  METHOD  OF  CHARACTERISTICS 


form  of  the  characteristic  equations  corresponding  to  characteristic 
lines  l”,  tl”,  III~  on  the  grid  In  the  left  side  of  Figure  22(b) 2 


I : 

‘^A 

- 

+ 

pc^A(v^  - 

il": 

Ma 

- 

plCb(a)A  - 

‘^B> 

IIl‘: 

^A 

- 

+ 

pAc^(w^  - 

(Aa) 

•c,  At  (Ab) 

hk  + 


At  (Ac) 


The  equations  for  axial 1 transverse » and  rotational  motion  of  the  front 
rigid  mass  are: 


dv 


®’a  " "f  dt 


f 


(5a) 


" ^A 


■Ma  + ^A^„ 
A A o 


dWj 

'f  "dT 

iiLi 

f dt 


(5b) 


(5c) 


where  the  subscript  f indicates  values  at  the  front  mass,  is  the 

front  mass,  is  axial  load,  F^  Is  lateral  load,  ^ Is  indicated 

In  Figure  21(a),  and  is  the  moment  of  inertia  of  the  front  mass  about 

the  point  on  the  axis  through  which  acts. 
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In  finite  difference  form  these  equations  are 


^ + "a  ■ ”f 


+ Oa  ■ "f 


/”a  - "p 

\ 2At 


•«A  + “ ^f 


(^k  ~ ‘^d\ 
\ 2At  / 


Because  the  same  values  of  v*  w,  and  u)  are  used  In  equations  (4)  and  (6)i 
compatibility  of  displacement  and  rotation  between  the  front  rigid  mass  and 
the  beam  Is  ensured, 

The  couples  motion  of  the  aft  mass  and  the  aft  end  of  the  beam  is 
governed  by  a similar  set  of  equations.  The  finite  difference  form  of 
the  characterkstlc  equations  for  the  aft  mass  Is 


I i Nj,  ~ v^)  - 0 


II  I + plc^  ((i),  - (Dy)  - C 


b ' Z X'  1) 


+ I 7 + 1 

III  : Q2  " Qy  “ “ -k"ACJ  \ At  (7c) 
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The  finite  difference  form  of  the  equations  of  motion  of  the  aft 


mass  are: 


-N  ■ m 
z a 


-Q  ■ tn 
z a 


M + Q Ji.  - I 
z ^z  2 a 


\ 2At  ) 

~ \ 

\ 2At  / 

\ 2At  I 


where  the  subscript  a refers  to  values  at  the  aft  massi  m is  the 

z 

aft  mass,  I is  indicated  in  Figure  21(b),  and  I is  the  moment  of 
2 ^ 

inertia  of  the  aft  mass  about  its  center. 

Each  set  of  six  simultaneous  linear  algebraic  equations,  (4)  with 
(6)  and  (7)  with  (8),  can  then  be  solved  for  their  six  unknowns  at 
points  A and  Z.  By  representing  values  of  the  variables  at  points 
C and  Y in  terms  of  their  values  at  points  A,  B,  and  D,  and  Z,  X, 
and  W,  the  solutions  can  also  be  expressed  in  terms  of  their  values  at 
the  corners  of  the  Cartesian  grid  of  Figure  22(b)  . The  complete  solution 
is  given  in  Appendix  B, 

Payload  Equations 

When  a penetrator  impacts  a target,  the  inertia  forces  from  an 
interior  payload  can  affect  the  response  of  the  penetrator  casing.  The 
effect  of  the  axial  inertia  can  be  eliminated  during  the  early-time 
response  simply  by  attaching  the  payload  with  a compliant  elastic  mount 
or  by  a weak  plastic  mount  such  as  crushable  foam.  The  transverse  or 
whipping  inertia  cannot  be  treated  in  a similar  way  because,  to  allow 
maximum  penetration,  the  outside  diameter  of  the  penetrator,  and  hence 
the  tolerance  between  the  penetrator  casing  and  the  payload,  must  be 
minimized.  Thus,  the  effect  of  the  payload  transverse  inertia  must  be 
examined . 





The  bending  ntlffnoss  of  the  payload  is  generally  much  less  than 
that  of  the  penetrator  casing.  Thus,  the  payload  is  modeled  as  a 
distributed  mass  with  no  bending  stiffness.  This  effect  was  modeled 
mathematically  by  increasing  the  density  from  p to  p'  in  the  equation 
of  translational  motion  of  the  beam.  The  corresponding  shear  velocity 

decreases  to  c'.  In  the  region  of  the  payload,  the  following  character- 
istic equation  must  be  used; 

nr-!  d[Q  + p"  A c'  w]  - -k  A G U)  dt  (9) 

Note  that  this  changes  the  location  of  points  H and  J in  Figure  22(b). 
Nevertheless,  these  equations  are  solved  in  the  same  way  as  the  beam 
equations  without  the  payload  but  with  appropriate  locations  for  H and 

J. 


RESPONSE  PREDICTED  IN'  SIMULATOR  EXPERIMENTS 

The  analysis  was  used  first  to  calculate  the  penetrator  strains 
measured  In  two  of  the  simulator  experiments  described  In  Section  3. 

The  structures  analyzed  were  a thick-walled  penetrator  model  (Test  64) 
and  a thln-walled  penetrator  model  (Test  66),  both  with  hemispherical 
noses;  the  dimensions  of  these  structures  are  given  In  Table  1.  The 
loads  measured  In  the  experiments  were  used  as  Input  for  the  analyses. 

To  account  for  the  uncertainty  in  the  point  of  application  of  the  re- 
sultant loading  force  In  the  simulator  experiments,  upper  and  lower 
bounds  on  the  response  were  calculated  for  extreme  locations  of  the 
point  of  load  application  for  each  experiment.  These  -omparlsons  not 
only  provided  a check  on  the  analysis  but  also  aided  In  the  understanding 
of  the  simulator  loading  technique. 

To  describe  the  loading  of  the  penetrator  In  a simulator  experiment, 

we  must  know  the  vertical  force  resultant  (l.e.,  the  force  parallel 

to  the  loading  piston  axis),  the  horizontal  force  resultant  F , the  angle 

of  inclination  9 of  the  penetrator  axis  relative  to  the  vertical  loading 

piston  axis,  and  the  points  of  application  of  F and  F (about  which 

V H 
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the  distributed  vertical  and  horizontal  forces  produce  no  moments) . 
From  the  Inset  In  Figure  23,  the  axial  and  transverse  force  resultants 
are,  respectively, 


Fa  - Fy  cos6  + Fy  sine 
and 

We  now  obtain  each  of  these  loading  parameters  from  the  experimental 
measurements.  The  vertical  force  is  obtained  from  the  load  cell  readings 
taken  during  the  test.  ‘ However,  the  load  cell  record  is  the  force  on 
the  bottom  side  of  the  sliding  disc,  which  differs  from  the  force  on  the 
pcnetrator  by  the  acceleration  force  of  the  disc.  The  strain  gage  records 
indicate  that  no  noticeable  force  is  applied  to  the  penetrator  during 
approximately  the  first  175  ysec  of  the  load  cell  record  [e.g.,  see 
Figure  18(a)]  when  all  the  gaps  In  Che  load  path  between  the  load  cell 
and  penetrator  are  closed  and  the  sliding  disc  begins  to  deform  under 
a small  force.  Therefore,  as  input  to  the  analysis,  we  omit  this  early 
portion  of  the  load  cell  record.  In  'Kkliclon,  subsequent  small  oscil- 
lations appear  in  the  load  cell  record.  These  are  caused  by  a vibration 
in  Che  piston.  This  vibration  is  not  transmitted  to  Che  penetrator  model 
througli  the  plastically  deforming  eliding  disc  as  evidenced  by  the  smooth 
strain  gage  records  [Figure  18(b) j.  Therefore,  the  oscillations  were 
also  smoothed  out  for  the  input.  Thus,  the  F^  curves  used  in 

the  analysis  of  these  two  tests  were  taken  as  shown  In  Figures  23  and  24, 
where  they  are  compared  with  the  actual  load  cell  records . 

A reasonable  assumption  concerning  F^^  is  that  the  ratio  of  F^ 
to  Fy  in  a given  experiment  is  constant  for  the  duration  of  the  loads. 
This  ratio  was  determined  from  high  speed  movies  of  the  penetrator  in 
free  flight,  during  which  the  center  of  gravity  of  the  penetrator  travels 
in  the  direction  of  the  resultant  impulse  (see  Section  3) . The  ratio 
Fjj/Fy  is  equal  to  the  slope  of  the  trajectory  of  the  center  of  mass  of 
the  penetrator,  which  could  be  measured  to  within  +10  percent.  For 
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Tests  64  and  66,  in  which  6 la  In  the  neighborhood  6f  20  degrees,  a 
10  percent  variation  in  Fj^  leads  to  only  a 2,5  percent  variation  in 
and  essentially  no  variation  in  F^,  In  Test  64,  ■ 0.141, 

and  in  Test  66,  “ 0.074. 

The  angle  of  inclination  of  the  penetrator  model  was  measured 
before  the  experiment  to  an  accuracy  of  ±15'.  In  Test  64  9 ■ 19'’30', 

and  in  Test  66  9 ■■  18^15'.  For  these  angles,  ±15'  variation  leads 

to  about  ±1  percent  variation  in  and  essentially  no  variation  in 
F^.  Since  the  exact  point  of  application  of  the  load  is  unknown,  the 
response  was  predicted  for  two  bounds  on  this  location.  One  bound  on 
the  location  of  the  point  of  load  application  is  through  the  nose  tip 
of  the  penetrator  at  the  axis.  This  is  a left  bound  (for  the  configuration 
shown  in  Figure  8)  because  the  penetrator  models  were  inclined  to  the 
vertical  so  that  the  point  of  load  application  is  actually  to  the  right 
of  the  penetrator  axis.  The  bound  on  location  to  the  right  was  taken 
as  the  center  of  the  dent  in  the  sliding  disk  because  the  resulting 
motion  of  the  penetrator  was  to  the  right  (Figure  8) . 

For  Test  64,  Figure  25  shows  the  upper  and  lower  bounds  on  pre- 
dicted strains  along  with  the  strain  measured  at  strain  gage  8.  The  error 
bands  on  each  of  the  predicted  strains  correspond  to  uncertainty  in  F^^ 
and  9.  The  total  strain  calculated  for  the  two  extreme  points  of  load 
application  bound  the  initial  peak  measured  in  the  experiment.  Since  the 
edge  of  the  dent  was  so  near  the  axis  tip  (0.06  inch),  it  is  reasonable 
that  the  upper  bound  calculation  lies  closer  to  the  measured  response. 

It  also  indicates  that  the  predicted  strain  depends  significantly  on  the 
point  of  load  application.  iTie  bounds  on  the  actual  point  of  application 
of  the  resultant  load  In  this  experiment  could  be  decreased,  for  example, 
by  taking  the  left  bound  (Figure  8)  as  the  left  edge  of  the  dent  in  the 
sliding  disk. 

For  Test  66,  Figure  26  shows  the  upper  und  lower  bound  calculated 
strains  along  with  the  strain  measured  at  gage  5.  Again,  the  error 
hands  represent  variations  In  and  9.  In  this  case,  the  first  peak 

in  strain  is  predicted  very  closely  by  the  lower  boucd  calculation, 
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FIGURE  2B  PREDICTED  STRAINS  AND  MEASURED  STRAIN  IN  SIMULATOR 
TEST  b4  (THICK-WALLED  STRUCTURE) 
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TOTAL  STRAIN  AT  REAR  LOCATION  - percent 


Indicating  that  was  applied  near  the  center  of  the  dent . This  is 
reasonable  since  the  dent  in  this  test  was  relatively  small  (0.3-lnch 
diameter) . 

These  comparisons  show  that,  with  accurate  knowledge  of  the  loads 
on  the  penetrator,  the  analysis  Is  capable  of  predicting  penetrator 
response  with  reasonable  accuracy.  It  also  shows,  however,  that  con- 
siderable differences  arise  from  uncertainty  In  the  point  of  application 
of  the  loads.  Thus,  although  the  predicted  response  is  accurate  for 
a given  input,  the  correlation  with  experiments  is  limited  by  the  quality 
of  the  experimental  load  data. 

APPLICATION  OF  THE  ANALYSIS 

Using  the  analysis,  we  now  determine  the  effects  of  eccentricity 
factor,  pulse  shape,  rise  time,  end  masses,  and  lateral  payload  on  pene- 
trator response  and  nose  rotation  as  it  affects  terradynamics. 

Eccentricity  Factor 

As  further  verification  of  the  applicability  of  the  analysis  and 
to  estimate  the  ratio  of  lateral  lo  axial  load  (called  the  eccentricity 
factor  n)  developed  on  Impact  at  an  angle  of  attack,  we  applied  the 
analysis  to  the  reverse  ballistics  Test  D-1  performed  at  AVCO.^  The 
eccentricity  factor  was  determined  by  adjusting  the  axial  and  lateral 
load  amplitudes  so  that  the  predicted  strain  response  best  matched  the 
measured  strain  response. 

The  penetrator  structure  used  in  AVCO's  Test  D-1  la  shown  in 
Figure  27.  This  structure  was  Impacted  by  a 15 . 2-lnch-diameter 
(38.6-cm-dlameter)  mortar  projectile  at  1500  ft/sec  (457  m/sec)  at 
a 5 degree  angle  of  attack.  The  strain  data  used  here  were  taken  from 
the  three  stations  Indicated  in  Figure  27.  From  these  strain  measure- 
ments, we  computed  the  axial  and  bending  strains  at  these  three  locations 
on  the  structure  at  each  of  three  times.  The  nine  axial  and  nine  bending 
strain  data  points  are  shown  In  Figures  28(b} , 28(c),  and  28(d). 


The  Idealized  structure  used  to  calculate  strains  is  also  shown  in 
Figure  27.  The  nose  length  is  4.5  inches  (11,4  cm).  In  the  structure 
tested,  this  length  is  the  distance  between  the  point  of  the  nose  and 
the  station  where  the  diameter  reaches  77  percent  of  the  diameter  of  the 
aft  cylindrical  section.  The  shapes  of  both  axial  and  lateral  load-time 
histories  were  taken  as  bilinear,  with  the  rise  time  of  250  ysec  cor- 
responding to  the  4.5-inch  nose  length  divided  by  the  Impact  velocity 
of  1500  ft/sec.  In  the  analysis,  the  loads  were  applied  at  the  center 
of  mass  of  the  rigid  nose  of  the  Idealized  structure. 

The  axial  and  lateral  load  amplitudes  that  give  best  agreement  be- 
tween measured  and  calculated  strains  were  determined  independently. 

For  example,  the  axial  response  was  first  calculated  for  an  axial  load 
of  unit  amplitude.  Then,  for  each  axial  strain  data  point,  the  unit 
load  was  weighted  such  that  the  calculated  axial  strain  agreed  with  the 
measured  strain.  The  axial  load  amplitude  that  best  matched  the  axial 
strain  data  was  taken  as  the  average  of  the  nine  such  weighting  factors; 
this  gives  ■ 195,000  pounds.  By  the  same  method,  we  found  that  a 
lateral  load  amplitude  of  53,000  pounds  best  matched  the  bending  strain 
data. 

These  load  histories  and  resulting  strain  distributions  are  shown 
in  Figure  28.  In  general,  good  agreement  was  obtained  between  the  calcu- 
lated and  measured  strains.  It  is  concluded  that  this  analysis  is 
capable  of  predicting  the  shape  of  the  strain  distribution  and  its  time 
history  in  a penetrator  under  angle-of-attack  Impacts.  For  this  5 degree 
angle  of  attack,  the  eccentricity  factor  Is  n - 53,000/195,000  ■ 0.27. 
Thus,  the  eccentricity  angle,  tan'^n  “ 15  degrees,  is  greater  than  the 
angle  of  attack  a by  a factor  of  about  3.  This  result  will  be  used 
in  Section  5 in  the  calculation  of  critical  Impact  curves. 

The  effect  of  varying  the  eccentricity  factor  was  studied  by  applying 
loads  of  different  eccentricity  factors  to  another  structure.  Since  only 
relative  changes  in  the  response  are  of  Interest,  the  structure  analyzed 
had  only  the  gross  characteristics  of  penetrators.  Specifically,  we 
analyzed  a simple  steel  tube  with  a ratio  of  length  to  outside  diameter 
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of  A/d  ■ 6 and  a ratio  of  radius-to-thlckness  of  a/h  ■ 4.  The  axial 
load  history  is  shown  in  Figure  29(a).  Calculations  were  made  for  seven 
different  eccentricity  factors  ranging  in  value  from  0 to  0.3. 

The  stress  histories  at  Station  x **  0.4828  for  the  tensile  surface 
(Side  A)  and  the  compressive  svirface  (Side  B)  of  the  structure  for  three 
different  load  ratios  are  shown  in  Figure  29(b).  Table  2 presents  the 
peak  tensile  and  compressive  stresses  and  the  times  and  locations  at 
which  they  occur.  The  time  and  location  of  peak  stress  change  mono- 
tonlcally,  but  only  slightly,  over  the  range  of  load  ratios,  resulting 
in  a nearly  linear  variation  in  peak  tensile  and  compressive  stresses 
with  load  ratio,  as  shown  in  Figure  29(c).  Thus,  for  a particular 
loading  rate  and  length  of  structure,  the  time  and  location  of  peak 
stresses  do  not  change  significantly,  and  their  magnitude  varies 
directly  with  load  ratio  for  ratios  of  10  percent  to  30  percent, 

Table  2 

DEPENDENCE  OF  MAXIMUM  STRESS  ON  TRANSVERSK/AXIAL 
FORCE  RATIO  (A/d  - 6,  a/h  " 4) 


Fl/Fa 

0 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

Maximum  tensile 

0 

0 

0.237 

0.542 

0.850 

1.173 

1.479 

stress/o 

8 

Location  - x 

— 

0.52 

0.48 

0.48 

0.45 

0.45 

Time  - t 

— 

“ 

4.41 

4.50 

4.50 

4.50 

4.50 

Maximum  compressive 
stress/a 

s 

-1 

-1.099 

-1.356 

-1.621 

-1.902 

-2.203 

-2.479 

Location  - x 

0 

0.21 

0.24 

0.31 

0.31 

0.34 

0.34 

Time  - T 

3.33 

3.75 

3.75 

3.75 

3.75 

3.84 

3.84 

Pulse  Shape 


The  effect  of  the  shape  of  the  loading-time  curve  on  penetrator 
response  was  studied  by  applying  loads  of  three  different  time  histories 
to  the  same  structure.  The  axial  load  time  histories  are  shown  In 
Figure  30(a);  in  each  case  the  ratio  of  the  lateral  load  to  the  axial 
load  was  0.2,  The  structure  was  the  same  simple  steel  tube  used  to 
study  the  eccentricity  factor  (A/d  ■ 6,  a/h  ■ 4). 

For  the  three  load  histories  the  corresponding  axial  and  (tensile) 
bending  stress  histories  at  Station  x ■ 0,2759  are  shown  in  Figure  30(b). 

At  T ■ 3r  that  is,  after  three  transits  of  a longitudinal  wave  over  the 
length  of  the  structure,  Figure  30(a)  shows  that  loading  Force  A has 
Just  reached  its  peak  amplitude,  Force  B Is  about  14  percent  less  than 
Force  A,  and  Force  C is  about  24  percent  less  than  Force  A,  At  r 3.5759  , 
the  time  required  for  the  disturbance  from  the  load  at  t ■ 3 to  arrive 
at  the  X " 0.2759  station.  Figure  30(b)  shows  that  the  bending  strei^s 
for  Load  B Is  about  10  percent  less  than  bending  stress  for  Load  A and 
that  bending  stress  for  Load  C is  about  20  percent  less  than  bending 
stress  for  Load  A.  Axial  stress  for  Load  B is  about  19  percent  less 
than  axial  stress  for  Load  A,  and  axial  stress  for  Load  C is  about  25 
percent  less  than  axial  stress  for  A.  Thus,  for  these  loading  curves, 
having  rise  times  of  three  to  five  transit  times,  the  stresses  generally 
follow  the  loading,  and  the  effect  of  the  detailed  shape  of  the  loading 
curve  on  the  stresses  is  small. 

Rise  Time 

The  effect  of  loading  rise  time  on  penetrator  response  was  studied 

by  applying  loads  with  the  different  rise  times  as  shown  in  Figure  31(a) 

and  with  an  eccentricity  factor  of  0.2.  The  structure  was  the  simple 

steel  tube  analyzed  previously.  The  bending  and  axial  stress  histories 

at  Station  x “ 0.4828,  normalized  with  respect  to  o ■ F /A,  corresponding 

B A 

to  each  loading  curve  are  Siiown  in  Figure  31(b). 
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(a)  LOAD  HISTORIES  (rj  » 0.20) 


r - C|,t/K 

(b)  AXIAL  AND  BENDING  STRESS  HISTORIES  AT  x “ x/K  =■  0.2759 

MA-3091-112A 

FIGURE  30  RESPONSE  OF  A SIMPLE  STEEL  TUBE  TO  DIFFERENT  LOAD 
HISTORIES  (£/d  “ 6,  a/h  - 4) 
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The  largest  bending  stress  value  (o./a  ■ 1.51)  Is  produced  by  the  ; 

D S 

loading  of  zero  rise  time  (step  Input) . In  order  of  Increasing  loading  ' 

rise  times,  the  peak  bendlt.g  stresses  o./a  produced  are  1.48,  1.40, 

1.31,  and  1,20.  As  the  rise  time  Increases  further,  the  peak  stress  at 
this  location  approaches  the  steady-state  value  of  ■ 0.793. 

Furthermore,  the  time  at  which  the  peak  occurs  Increases  with  rise  time. 

For  the  axial  stress,  the  peak  stress  at  a location  is  determined 

by  whether  the  loading  rise  time  is  less  than  or  greater  than  the  time 

required  for  a reflected  tensile  stress  wave  to  arrive  at  that  location 

from  the  aft  free  end.  Fur  the  loading  curves  with  rise  times  of  less 

than  one  transit  time,  the  maximum  value  of  a /o  is  -1.00.  For  the 

a s 

loading  curves  wl^-h  rise  times  of  more  than  one  transit  time,  the  maximum 
stress  is  o /o  » -0.617.'  The  steady-state  value  at  this  location  is 

‘j'i. 8 S 

0 /a  ■ -0.517. 
a 8 

Thus,  we  conclude  that  the  largest  effects  of  increasing  the 
loading  rise  time  are  to  increase  the  times  at  which  peak  bending  occurs 
and  to  vary  the  axial  stress  ‘'or  loading  rise  times  near  one  transit 
time.  For  rise  times  larger  than  two  transit  times,  the  variation  in 
loading  rise  time  han  less  effect  on  the  peak  bending  stress  and  the 
peak  axial  stress. 

Enrl  Masses 

The  effect  of  including  the  forward  and  aft  rigid  end  masses  in 
the  calculation  of  penetrator  response  was  determined  by  analyzing  the 
structure  shown  in  Figure  32(a).  The  central  portion  nf  the  structure 
is  the  steel  tube  analyzed  in  previous  calculations.  The  front  conicul 
mass  has  a length  equal  to  the  outside  diameter  d,  and  the  aft  cylin- 
drical mass  has  a length  equal  to  1/3  the  diameter.  Thus,  the  forward 
and  aft  masses  are  equal,  each  being  about  12.7  percent  of  the  mass  of 

the  hollow  uection  Four  versions  ot  this  structure  were  analyzed: 

without  end  masses,  with  a front  mass,  with  an  aft  mass,  and  with  both 
front  and  aft  masses.  In  each  case  the  eccentricity  factor  Is  ti  - 0,2 

and  the  rise  time  Is  i ■ 2 . 5 as  shown  in  the  top  of  Figure  32 (b).  For 
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a structure  with  a front  massi  the  lateral  load  Is  applied  through 

the  center  of  gravity  of  the  front  mass.  For  a structure  without  a 
front  mass,  the  lateral  load  is  applied  at  the  front  end  of  the  hollow 
section. 

The  predicted  bending  and  axial  stress  histories  at  Station  x " 
yi/i  • 0.4138  are  shown  in  Figure  32(a)  and  (b).  This  station  is  the 
location  of  maximum  bending  stress  for  the  structure  without  end  masses. 

The  locations  of  maximum  bending  stress  for  the  other  cases  are  x " 0.3793 
with  the  front  mass  x ■ 0.4483  with  aft  mass,  and  x ■ 0.4138  with 
both  masses.  These  are  adjacent  locations  on  the  finite  difference  grid 
for  the  structure.  Thus,  the  presence  of  end  masses  has  virtually  no 
effect  on  the  location  of  the  maximum  bending  stress. 

The  magnitude  of  the  peak  bending  stress,  however.  Is  affected  by 
the  end  masses.  There  is  a 31  percent  reduction  when  the  front  mass 
is  included,  a 12  percent  increase  when  the  aft  mass  is  Included,  and  a 
19  percent  reduction  when  both  masses  are  Included.  Including  either  the 
front  or  aft  mass  increases  the  time  at  which  the  peak  response  occurs 
by  about  14  percent.  Including  both  masses  increases  this  time  by  about 
28  percent. 

Similarly,  for  the  axial  stress,  there  is  a 14  percent  reduction  in 
peak  stress  when  the  front  mass  is  included,  an  insignificant  change  in 
peak  stress  when  the  aft  mass  is  included,  and  about  a 14  percent  re- 
duction in  peak  stress  when  both  masses  are  included.  A reduction  of 
about  14  percent  in  frequency  of  response  occurs  with  inclusion  of  the 
front  or  aft  mass;  Inclusion  of  both  masses  practically  eliminates  the 
stress  oscillations. 

Thus,  we  conclude  that  inclusion  of  the  front  mass  reduces  the  peak 
axial  and  bending  stresses,  inclusion  of  the  aft  mass  increases  the 
bending  stress,  and  inclusion  of  both  masses  produces  stresses  interme- 
diate between  these  extreme  values.  Inclusion  of  either  mass  increases 
the  time  at  which  the  peak  stress  occurs. 


Lateral  Payload  Inertia 


The  effect  of  Including  the  lateral  payload  inertia  waa  determined 
by  analyzing  the  structure  shown  in  the  inset  in  Figure  33  with  and  with- 
out the  payload.  The  cylindrical  section  of  the  casing  is  made  of  steel 
of  density  0.281  Ib/in^  (7773  kg/m^)  and  is  geometrically  similar  to 
that  of  a shallow  penetrator.  The  payload  has  a density  of  0,062  Ib/ln^ 
(1713  kg/tn^)  and  a mass  that  is  78  percent  of  the  mass  of  the  cylindrical 
section.  The  mass  of  the  front  nose  is  38  percent  of  the  mass  of  the 
cylindrical  section.  Since  the  wall  of  this  structure  is  as  chin  as  a 
penetrator  casing  is  likely  Co  be  made»  this  configuration  emphasizes 
the  effect  of  the  payload  lateial  inertia. 

The  calculated  bending  stress  histories  shown  in  Figure  33  are 
similar  with  and  without  the  payload.  The  maximum  difference  is  about 
14  percent  and  occurs  at  the  time  of  the  maximum  bending  strain. 

Further  response  calculations  do  not  include  the  effect  of  lateral 
payload  inertia  for  two  reasons.  Firsti  the  effect  on  total  strain  is 
small.  In  an  actual  penetrator  the  payload  could  be  made  to  slide 
relative  to  the  casingi  allowing  the  nose  to  decelerate  the  payload  so 
that  the  axial  component  of  stress  in  the  casing  would  not  be  altered  by 
the  payload.  Therefore,  the  effect  of  the  payload  on  total  (bending 
plus  axial)  stress  is  even  less  (by  about  a factor  of  2)  than  the  effect 
on  bending  stress.  Second,  the  critical  impact  curves  of  Section  5 are 
most  useful  for  comparing  differences  in  peak  response  stress  among 
several  structures.  The  difference  in  peak  response  stress  from  one 
structure  to  another  is  similar  with  or  without  the  payload. 

Nose  Rotation 

We  also  determined  the  rotation  of  the  nose  to  assesB  whether  de- 
formation of  the  structure  could  alter  the  terradynamlcs  and  therefore 
the  load.  For  the  simple  tube  structure  without  end  masses  [Figure  31(a)], 
the  rotation  of  the  front  of  the  structure,  relative  to  Its  orientation 
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FIGURE  33  BENDING  STRAIN  AT  x - 0,27  VERSUS  TIME  FOR  A 

PENETRATOR  WITH  AND  WITHOUT  PAYLOAD  (C/d  » 3.6, 
a/h  * 10.0,  pVp  - 1.9B5,  t,  - 6.5,  r?  - 0.2) 
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when  the  load  was  applied,  Is  only  0.032  degree  at  the  time  peak  bending 
occurs.  Thus,  It  Is  concluded  that  during  penetration  this  rotation  is 
not  great  enough  to  affect  the  loading  on  the  nose. 


CONCLUSIONS 

Comparison  of  calculated  response  with  that  measured  In  the  simulator 
experiments  and  In  a reverse  ballistics  test  showed  that  this  simple 
beam-mass  theory  Is  sufficient  for  predicting  the  strain  response  In 
penetrator  structures.  The  parameters  that  have  the  largest  effect  on 
the  peak  compressive  strain  are  the  eccentricity  factor,  the  rise  time, 
and  the  end  masses. 


I 
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5 CRITICAL  IMPACT  CURVES 
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The  analysis  was  used  to  calculate  the  response  of  penetrator 
structures  ranging  from  designs  for  deep  penetration  (large  i/d, 
small  a/h)  to  those  for  shallow  penetration  (small  i/d,  large  a/h) . 
Calculations  are  also  made  for  penetrators  with  Intermediate  values  of 
i/d  and  a/h. 

For  design  purposes > the  most  useful  Information  la  the  relationship 
between  the  Impact  parameters  and  the  response  parameters . Then  design 
tradeoffs  can  be  made  between  the  penetrator  structure  and  the  Impact 
conditions.  For  example)  for  a given  penetrator  (l.e.i  a given  allowable 
peak  response  stress)  and  a given  target,  a tradeoff  can  be  made  between 
Impact  velocity  and  angle  of  attack. 

The  complete  impact-response  relationship,  including  terradynamlcs. 

Is  complex  and  not  well  understood  for  many  Impact  conditions.  However, 
the  procedure  for  determining  this  relationship  can  be  simplified  and  some 
Insight  Into  the  Impact-response  relationship  can  be  gained  by  Introducing 
an  appropriate  set  of  parameters  that  describe  the  loading  on  the  pene- 
trator. The  Impact-load  relationship  (from  terradynamlcs)  and  the  load- 
response  relationship  (from  structural  response)  can  be  treated  separately 
and  then  combined  to  produce  the  impact-response  relationship  sought. 
Because  terradynamlcs  Is  not  well  understood  we  must  make  reasonable 
assumptions,  based  on  currently  available  Information,  about  the  Impact- 
load relationship.  The  load-response  relationship  can  be  determined  from 
the  structural  analysis  described  In  Section  4. 

CRITICAL  IMPACT  CURVES  FOR  SIMPLE  TUBE  STRUCTURES 

This  procedure  is  Illustrated  below  for  the  simple  tube  structure 
with  i/d  ■ 6 and  a/h  ■ 4 analyzed  in  Section  A.  The  load-response 
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relationship  Is  discussed  first,  then  the  Impact-load  relationship. 

These  two  relationships  are  combined  to  produce  the  Imp act -response 
relationship. 

Load-Response  Relationship 

The  dominant  loading  parameters  are  the  rise  time  the  peak 

axial  force.)*'  F,  and  the  ratio  n of  lateral  load  to  axial  load.  The 

dominant  response  parameter  Is  the  peak  compressive  stress 

relationship  among  these  parameters  can  be  determined  from  the  results 

of  the  calculations  of  the  effects  of  rise  time  discussed  in  Section  3. 

Ihe  results  of  some  of  these  calculations  were  given  In  Figure  31,  which 

shows  the  bending  and  axial  stress  at  x ••  0.4828  for  n ” 0.2  and  for 

dimensionless  rise  times  2.55,  and  3.4. 

We  first  determine  the  normalized  peak  response  stress  as  follows.  For 

each  rise  time  and  value  of  n we  find  a , the  total  compressive 

max 

stress  that  Is  maximum  In  space  and  time  during  the  time  of  Interest, 
and  then  normalize  it  with  respect  to  the  amplitude  of  the  normal  stress 
a on  the  loaded  end.  For  example,  In  Figure  31(b)  for  t ■ 0 and 

6 i 

n ■ 0.2,  the  maximum  value  of  a,  /a  is  1.51  at  t ■ 3 and  x “ 0.4828 

D 8 

and  the  maximum  value  of  Is  1.0  over  a range  of  values  of  x* 

Thus,  the  maximum  value  a h Is  2.51  at  t ■ 3 and  x ■ 0.4828. 

max  8 

Curves  similar  to  those  shown  In  Figure  31  can  also  be  used  to 

determine  c /a  for  different  values  of  n.  Since  the  axial  and 
max  8 

bending  responses  are  independent,  the  total  stress  for  different  values 
of  n can  be  computed  by  weighting  the  bending  component  by  the  appro- 
priate factor.  For  example,  for  n “ 0.1,  the  bending  stress  shown  in 
Figure  31  (for  which  n ■ 0.2)  is  multiplied  by  0.1/0. 2 ■ 0.5  before 
being  added  to  the  axial  stress. 


For  the  small  angles  of  attack  considered  here,  the  amplitude  of  the 
axial  load  is  approximated  by  the  amplitude  of  the  load  F parallel 
to  the  penetrator  velocity  vector. 
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The  normalized  peak  responee  stress  as  a function  of  rise  time  for 
several  values  of  n are  plotted  In  Figure  34.  Note  that  this  load- 
response  relationship  applies  for  any  Impact  conditions  that  produce 
the  appropriate  values  of  T^,  F,  and  n* 

Impact-Load  Relationship 

The  dominant  Impact  parameters  are  the  Impact  velocity  V,  the  angle 
of  attack  a I and  the  target  impedance*''  Z.  Thus,  the  Impact-load 
relationship  relates  the  loading  parameters  F,  and  n to  V,  ot, 
and  Z.  As  discussed  previously,  the  first  loading  parameter,  the  rise 
time  T , is  well  approximated  by  the  submersion  time  of  the  tapered 
penetrator  nose  length  SL^t 

% - c,  C9) 

r bn 

The  dependence  of  the  loading  force  F on  the  Impact  parameters  is 
not  as  well  understood.  However',  a reasonable  approximation  is  that  F 
Is  proportional  to  the  impact  velocity  V, 

F - b(Z)  V (10) 

where  the  proportionality  factor  b(Z)  depends  on  the  target  material. 

For  normal  Impacts  this  relationship  has  been  verified  both  experimentally 
and  theoretically  for  soil  targets;  It  also  appears  to  hold  for  rock 

t 

targets  but  fewer  experimental  data  are  available. 

The  dependence  of  the  load  ratio  n on  the  Impact  parameters  Is 
complex  and  has  not  been  adequately  Investigated.  However,  some  infor- 
mation is  available.  For  normal  impacts,  that  is,  for  a ■ 0,  tan“i  n ■ 0. 
For  angle-of-attack  Impacts,  additional  information  can  be  inferred  from 
the  AVCO  reverse  ballistics  test  data  and  structural  calculations  shown 


The  Impedance  Z Includes  any  target  properties  affecting  the  loading 
parameters,  including  moduli,  strength,  density,  and  viscosity  and 
frictional  coefficients. 

+ 

P.  F.  Hadala,  private  communication,  January  1977. 
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FIGURE  34  DYNAMIC  AMPLIFICATION  FACTOR  VERSUS  RISE  TIME  FOR 
A SIMPLE  STEEL  TUBE  <8/d  ■ 6.  a/h  - 4) 


In  Figure  28.  For  that  testt  a ■■  5 degrees  and  the  beat  agreement 
between  measured  and  predicted  strain  was  obtained  for  n " 0.267  or 
tan"^n  * 15  degrees.  Thus  for  this  case,  tan"^n  “ So.  We  now 
make  the  reasonable  assumption  that  the  eccentricity  angle  Increases 
linearly  with  angle  of  attack;  that  Is,  we  assume  that 


over  the  range  of 


an"^  ■ 3a 
of  Interest. 


Imnact-Responae  Relationship  and  Critical  Impact  Curves 


The  impact-response  relationship  Is  now  found  by  relating  the  three 
parameters  and  n In  the  load-response  relationship  of 

Figure  34  with  the  Impact  parameters  V and  o,  using  equations  (9), 
(10),  and  (11). 

First,  Is  related  to  impact  velocity  by  using  equation  (9). 

This  requires  specifying  the  nose  length  Here  we  have  chosen 


1.5d.  Then 


Second,  a ^ a Is  related  to  the  impact  force  F through 

max  8 

equation  (10).  However,  so  that  the  results  will  apply  to  any  targets 
for  which  a linear  force-velocity  relatio.iship  is  reasonable,  we  normalize 
the  peak  response  stress  with  respect  to  the  average  loading  pressure 
P ■ F/A"  at  a particular  velocity,  where  A"  is  the  frontal  area  of 
the  penetrator.  For  this  normalization  we  choose  P^OOO’  loading 

pressure  generated  at  an  Impact  velocity  of  1000  fc/sec.  The  peak 
response  stress  normalized  with  respect  to  P^^qqq  Is  related  to  the 
normalized  stress  of  Figure  34  through  the  identity 


1000 


F(IOOO) /A 
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Substitution  of  equation  (10)  Into  equation  (12)  and  rearrangement  gives 


0 

max 

^1000 


(14) 


where  V is  In  ft/sec.  For  the  simple  steel  tube,  k'/k  • 2.29. 

Finally,  the  n ■ constant  curves  of  Figure  34  are  relabeled  as 
ot  " constant  curves  according  to  equation  (11).  The  resulting  Impact- 
response  relationship  Is  shown  In  Figure  35.  Note  that  although  this 
Impact-response  relationship  relies  on  two  assumptions  concerning  the 
impact-load  relationship,  the  load-response  relationship  Involves  no  such 
assumptions.  Thus  when  more  Information  about  the  impact-load  relation- 
ship Is  available,  the  procedure  developed  here  can  be  repeated  with  the 
same  load-response  relationship. 


A more  useful  form  of  the  impact-response  relationship  Is  found  by 
constructing  cross  plots  from  Figure  35  for  which  the  normalized  peak 
response  stress  is  constant,  as  shown  In  Figure  36,  These  curves  are 
called  critical  Impact  curves  and  give  the  tradeoff  between  Impact 
velocity  and  angle  of  attack  at  the  maximum  capability  of  the  penetrator. 

That  is,  for  a given  penetrator  (fixed  fJjjjgjj)  and  target  (fixed  Pj^qqqi 
determined  experimentally  or  analytically),  the  curve  for  the  appropriate 
ratio  of  o gives  the  combination  of  allowable  values  of  V and 

a.  For  example,  for  '^niax^^lOOO  “ penetrator  is  to  survive 

at  angles  of  attack  as  large  as  5 degrees,  the  Impact  velocity  must  not 
exceed  1400  ft/sec. 


The  critical  Impact  curves  can  also  be  used  to  select  the  penetrator 
material,  on  the  basis  of  yield  stress,  for  a given  range  of  V and  a 
within  which  the  penetrator  must  function.  For  example,  If  for  a given 
system  and  target  the  maximum  Impact  velocity  is  2000  ft/sec  and  the 
maximum  angle  of  attack  Is  5 degrees,  a material  with  a yield  strength 
of  at  least  10  P^^qqq  Is  needed.  In  a similar  way  the  critical  Impact 
curves  of  Figure  36  could  be  used  to  select  targets  for  which  a given 
system  (l.e.,  specified  combinations  of  a^,  V,  and  a)  could  be  used. 
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Critical  Impact  curves  are  perhaps  most  useful  for  cooiparlng  the 
performance  of  different  penetrator  structures.  This  Is  Illustrated 
In  the  next  section. 

CRITICAL  IMPACT  CURVES  FOR  PENETRATOR  STRUCTURES 

Critical  impact  curves  were  developed  fur  the  four  penetrator 
structures  shown  In  Figure  37.  Structure  B Is  similar  to  that  already 
used  for  deep  penetration  into  hard  targets.  Structure  C Is  similar  to 
that  proposed  for  shallow  penetration.  Structures  A and  D were  analyzed 
to  determine  the  effects  of  changing  i/d  (along  the  rows  of  Figure  37) 
or  a/h  (along  the  columns  of  Figure  37). 

The  curves  of  the  normalized  peak  response  stress  versus  rise  time 
for  these  four  structures i with  n as  a parameter,  are  shown  In  Figures 
38  through  41.  Two  trends  are  apparent  in  these  curves.  First,  for 
two  structures  of  the  same  length  and  outside  diameter,  the  structure 
with  the  thicker  wall  has  a greater  spread  among  the  curves  (compare 
Figures  38  with  40  and  39  with  41) . This  is  because  material  added  to 
the  inside  wall  Increases  the  lateral  inertia  proportionately  more  than 
the  lateral  (bending)  stiffness.  Thus,  the  bending  stress,  and  therefore 
the  total  stress.  Is  greater  for  the  thicker  walled  structure.  Second, 
for  two  structures  of  the  same  thickness,  the  shorter  structure  has 
curves  with  a steeper  slope  for  small  t'  and  a change  In  slope  at 
t'SZO  (compare  Figure  38  with  39  and  40  with  41).  For  t"  > 20,  the 
response  of  the  shorter  structures  Is  quasl'-steady ; that  Is,  the  peak 
stress  depends  on  the  load  amplitude  but  not  the  rise  time.  In  contrast, 
for  the  rise  times  considered  here,  the  peak  stress  In  the  longer  struc- 
tures depends  on  both  rise  time  and  load  amplitude. 


Since  all  four  structures  have  the  same  diameter,  the  abscissa  coordinate 
■ c^^t^/d  along  the  top  of  each  graph  is  on  the  same  scale  In  all  four 

figures.  Thus,  we  can  make  direct  comparisons  of  the  dependence  of  the 
peak  response  on  rise  time  for  the  four  structures. 
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FIGURE  37  PENETRATOR  STRUCTURES  ANALYZED 
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FIGURE  38  DYNAMIC  AMPLIFICATION  FACTOR  VERSUS  RISE  TIME  FOR 
STRUCTURE  A (K/d  - 3.6,  a/h  = 2.  - 1.6d) 


FIGURE  40 
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The  critical  load  curves  for  these  four  structures  are  shown  In 
Figures  42  through  45.  The  trends  Identified  in  the  previous  curves 
of  peak  response  versus  rise  time  carry  over  to  this  set  of  curvcE. 

First,  for  the  thicker-walled  structures,  the  greater  spread  in  the 
curves  of  peak  response  versus  rise  time  give  steeper  critical  i- pc'L' 
curves  (compare  Figure  42  with  44  and  43  with  45).  Second,  fi  . the 
shorter  structures  and  for  small  t',  the  steeper  slopes  of  the  curves 
of  peak  response  stress  versus  rise  time  give  less  steep  critical  Impact 
curves  at  high  velocity  (compare  Figure  42  with  43  and  44  with  45). 

Critical  Impact  curves  similar  to  those  in  Figures  42  through  45  can 
be  used  to  plan  experiments.  For  example,  the  critical  load  curves 
for  structure  B (Figure  43)  shows  that,  for  impacts  at  about  1000  ft/sec 
and  a 2-degree  angle  of  attack,  the  severity  of  the  response  is  very 
sensitive  to  the  angle  of  attack.  Thus  it  may  be  Important  to  be  able 
to  control  the  angle  of  attack  in  the  experiment  or  at  least  be  able  to 
measure  it  accurately.  In  contrast,  Structure  C (Figure  44)  is  more 
sensitive  to  impact  velocity.  Thus  in  testing  Structure  C it  is  more 
important  to  control  or  measure  Impact  velocity  than  angle  of  attack. 

Critical  Impact  curves  can  also  be  used  to  interpret  results  of 
experiments  or  more  detailed  load  and  structural  response  calculations 
(e.g,,  a finite  element  code  prediction).  A particular  experiment  or 
detailed  calculation  gives  a single  point  on  a critical  Impact  curve. 
Although  such  a data  point  represents  a more  complete  response  descrip- 
tion than  for  the  beam-mass  model  used  here,  many  points  are  needed  to 
determine  the  shape  of  the  curve.  The  curves  calculated  w.1  i;h  the  beam- 
mass  model  can  be  used  to  determine  this  shape.  Thus  the  appropriate 
values  of  V and  a can  be  selected  to  minimize  the  data  points 


As  for  the  simple  steel  tube,  these  curves  were  generated  from  curves 
\aJ^1000  (not  shown). 
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FIGURE  44  IMPACT  VELOCITY  VERSUS  ANGLE  OF  ATTACK  FOR  CONSTANT 
PEAK  RESPONSE  STRESS  IN  STRUCTURE  C (e/d  = 3,6,  a/h  - 10.0) 
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needed  to  define  the  curve.  For  example,  for  Structure  B small  incre- 
ments in  a and  larger  increments  in  V should  be  made  to  efficiently 

define  the  o ■ 10  curve  for  small  a, 

max  1000 

The  critical  impact  curves  are  perhaps  most  useful  for  making  design 

tradeoffs  among  candidate  penetrator  structures.  For  example,  Figure  A6 

shows  the  o /‘P-mnn  ■ 10  curves  for  each  of  the  four  structures  analyzed. 

max  1000 

The  curves  for  Structure  B (deep  penetrator)  and  for  Structure  C (shallow 
i penetrator)  pass  through  a similar  region  in  the  V-a  plane  and  inter- 

sect at  V " 1500  ft/sec  and  a » 5.7  degrees.  However,  the  curve  for 
j Structure  B (deep  penetrator)  is  steeper  and,  for  small  angles  of  attack, 

! this  structure  can  withstand  greater  impact  velocities.  Both  these  pene- 

f 

I trators  can  be  made  stronger  by  increasing  only  the  wall  thickness 

i 

f (from  that  of  Structure  C to  that  of  Structure  a)  or  by  decreasing  the 

I length  (from  that  of  Structure  B to  that  of  Structure  A).  Note  that 

I either  of  these  changes  reduces  the  payload  volume.  Also,  for  small 

, angles  of  attack,  the  velocity  range  can  be  extended  more  by  Increasing 

wall  thickness  than  by  decreasing  length.  Similar  conclusions  can  be 
‘ drawn  by  comparing  curves  for  Structures  B and  C with  that  of  Structure  D. 

i 

CONCLUSIONS 

A method  has  been  developed  for  characterizing,  in  terms  of  critical 
Impact  curves,  the  ability  of  penetrator  structures  to  withstand  impacts. 

The  principal  advantage  of  this  procedure  is  that  it  characterizes  the 
I response  over  a range  of  impact  conditions.  This  allows  tradeoffs  to 

I be  made  among  structural  dimensions,  yield  strength,  and  target  character- 

I Istlcs  and  allows  selection  of  the  best  structure  from  a group  of  structures 

I for  a particular  application.  Critical  Impact  curves  can  also  be  used 

L 

I to  plan  and  Interpret  experiments  and  more  detailed  calculations. 
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FIGURE  46  IMPACT  VELOCITY  VERSUS  ANGLE  OF  ATTACK  FOR  CONSTANT 
PEAK  RESPONSE  STRESS  FOR  DIFFERENT  STEEL  PENETRATOR 
STRUCTURES  (ffmax^Pjoao 
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Appendix  A 
EXPERIMENTAL  DATA 


Table  A-1 

CALIBRATION  TESTS  WITHOUT  MODEL  PENETRATOR 


Test 

No. 

Date 

Charge 

Mass* 

(am) 

Initial 

Piston 

Dlsp. 

(cm) 

1 

3-12-76 

0.66*^ 

0.193 

2 

3-15-76 

0.66^ 

0.193 

3 

3-26-76 

0.543 

0.457 

4 

3-26-76 

0.543 

0.457 

5 

3-29-76 

1.086 

0.457 

6 

3-29-76 

1.086 

0.457 

7 

3-30-76 

1.086 

0.457 

8 

3-30-76 

2.00 

0.457 

9 

3-30-76 

2.00 

0.762 

10 

4-02-76 

1.00 

0.305 

11 

4-02-76 

1.00 

0.152 

Vent 

Holes 


5,  790 

6.070 
4,140 
3,650 
7,760 
7,760 

7.070 
14,130 
11,720 

7,450 

10,340 


Peak 

Acc. 

_e 

Rise 

Time'^ 

(nsec) 

1,700 

100 

1,050 

180 

850 

170 

2,000 

150 

2,400 

110 

1,800 

135 

4,000 

120 

3,600 

240 

2,000 

145 

2,650 
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^For  tests  3 through  11,  a small  (0.02  to  0.05  gram)  booster  charge 
(Du  Pont  Detasheet)  was  used;  the  mass  listed  here  includes  that  of 
90/10  PTEN/mlcrosphere  charge  only. 

^Average  of  the  two  pressure  gage  measurements. 

^Time  to  reach  75%  of  the  peak  acceleration. 

'^Includes  large  (approximately  0.18  gram)  booster  charge  (Du  Pont  Detasheet). 

^Bad  gage  record. 

^Vent  hole  configuration:  top  pair  open;  middle  pair  open;  bottom  pair 
closed . 

®Vent  hole  configuration;  top  pair  open;  middle  pair  closed;  bottom  pair 
closed . 


Table  A-2 

CALIBRATION  TESTS  WITH  LONG  ROD 


Test 

No. 

Data 

Device  Setup* 
Charge  Initial  Platon 
Maat^  Dlaplacement 

Iri) icffi) 

Peak 

Preasuce^ 

(kPa) 

Peak  Strain 
at  Aft 
Statlon<^ 

(*) 

Remarka 

12 

4-06-76 

1.00 

0.152 

Detonator  misfired. 

13 

4-07-76 

1.00 

0.152 

11,200 

0.021* 

14 

4-15-76 

1,00 

0.152 

— 

— 

Detonator  misfired. 

15 

4-15-76 

1.00 

0.152 

“ 

“ 

Detonator  misfired. 

16 

4-16-76 

1.00 

0.152 

11,000 

0.020 

17 

4-16-76 

1.00 

0.081 

— 

„ 

Detonator  misfired. 

18 

4-16-76 

1.00 

0.081 

12,400 

0.021 

19 

4-19-75 

1.00 

0.081 

13,800 

0.027 

20 

4-19-76 

1.00 

0.081 

13,800 

0.027 

21 

4-20-76 

3.00 

0,094 

43,400 

0.103^ 

22 

4-21-76 

3.00 

0.094 

— 

— 

Detonator  misfired. 

23 

4-21-76 

3.00 

0,094 

46,100 

0.090 

24 

5-13-76 

3.00 

0.094 

48,900 

0.089 

25 

S-13-76 

3.00 

0.094 

46,500 

0.087 

26 

5-14-76 

3.00 

0.305 

34,500 

0.053 

Special  and  condition  taated. 

27* 

.4-27-76 

1.00 

0.094 

12,200 

0.0208 

15.2-cin-rBdlue  loaded  end. 

28 

5-28-76 

1.00 

0.094 

10,700 

0.014 

Load  cell  tested. 

29 

6-23-76 

1.00 

0.152 

— 

— 

No  gage  record!. 

30 

6-23-76 

1.00 

0.152 

9,100 

0.016 

Load  call  teatid. 

*Addl Clonal  device  aeCup  parameCeia  fot  Toata  12  through  30; 

Vent  hole  configuration;  top  pair  open;  middle  pair  open;  bottom  pair  cloaed, 
eKcept  for  TaaC  27  where  all  three  paira  were  open. 

Orifice  Area;  1.4SS  cm^ 

For  Xeata  12  through  30,  a amall  (0.02  to  0.06  gram  booater  charge  (Du  font  Detaaheet) 
wai  uaed;  the  maaa  llaCod  hers  Includes  that  of  90/10  PETN/mlcrosphere  charge  only. 

'^Average  of  the  two  pressure  gage  measuromenta. 

'^Average  of  three  or  four  gagea  from  aft  station  (except  as  noted). 

^Gage  5 only. 

^Gage  1 only. 

^Average  of  Gages  1,  2,  and  3. 
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Table  A-3 


MODEL  STRUCTURE  TESTS 


Model 

Penetrator*' 


Initial  Angle 
[degree  (radii 


Peak  Load 
(Ml 


Strains  at 
Aft  Station  (%) 

Bending^  Peak  Axlal^ 


n S A 

31 

6-24-76 

SS 

Normal 

f 

90 

m m 

0.022® 

35 

7-20-76 

ss 

Normal 

89.0 

-- 

0.018 

36 

7-22-76 

SS 

Normal 

89.0 

0.018 

37 

7-22-76 

ss 

Normal 

77.8 

-- 

0.015 

32 

7-07-76 

ss 

4°  (0.07) 

73.4 

0.006 

0.016 

34 

7-09-76 

ss 

5‘  (0.09) 

71.2 

0.015 

0.015 

33 

7-08-76 

ss 

8“  (0.14) 

75.6 

0.018 

0.017 

38 

8-03-76 

TkW 

Normal 

89.0 

-- 

0,029 

39 

8-03-76 

TkW 

Normal 

66.7 

— 

..Vi 

40 

8-04-76 

TkW 

Normal 

73.4 

-- 

0.027 

i 

8-04-76 

TkW 

Normal 

75.6 

-- 

0.025 

43 

8-06-76 

TkW 

7°10'  (0.125) 

71.2 

0,016 

0.023 

42 

8-05-76 

TkW 

7-20'  (0.128) 

70^ 

0.020 

0.021 

i ^5 

8-17-76 

TkW 

a^SO'  (0.148) 

69.0 

0.022 

0.028 

I 44 

8-16-76 

TkW 

e^AS'  (0.153) 

62.3 

0.020 

0.020 

! 

8-19-76 

TAl 

Normal 

46.7 

-- 

0.030 

1 49 

8-20-76 

TAl 

S^lO'  (0.108) 

35.6 

..1 

46 

8-18-76 

TAl 

S-n'  (0.109) 

44,5 

0.005 

0,027 

48 

8-20-76 

TAl 

7°0'  (0.122) 

51.2 

0.008 

0.032 

*Devlce  setup  parameters  for  Tests  31  through  49: 

Vent  hole  configuration:  top  pair  open;  middle  pair  open;  bottom  pair  closed. 
Charge  mass;  2,00  gram  90/10  PETN/ralcrosphere  charge.  No  booster  charge  was 

used  except  In  Test  31  where  a 0,03-gram  patch  of  Ou  Pont  Detasheet 
was  used  In  addition  to  the  2.00-gram  90/10  charge. 

Initial  piston  displacement:  0,132  cm. 

Orifice  area:  1.455  cm^. 

^Abbreviations:  SS--Bolld  steel  model;  TkW-- thicked-wulled  model;  TA1-- tapered  alumi- 
num model, 

'^Peak  vertical  load  measured  with  load  cell. 

*^Bendlng  strain  at  time  of  peak  axial  strain  (average  magnitude  of  two  gages). 

Average  of  four  gages, 

f 

Determined  from  pressure  record. 

^Average  of  three  gages. 

^Teat  Incorrectly  set  up. 

^Bad  gage  records. 
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of  four  gages. 


Referring  to  Figure  21,  the  finite  difference  form  of  the  equations 
for  the  front  end  are 


'’a  ■ ■ Vj  . 0 

(1) 

^A  ■ S - “"b(  ^ 2 

(2) 

/ \ J/U).  + U)-\ 

%■%*  - "c)  • 2 )“ 

(3) 

' "ff  °) 

(4) 

/w.  - w \ 

■'l  + ’a  • ‘l{  2A,  ) 

(5) 

“a  ^A-^0  *■  ^f(  2At  ) 

(6) 

Let 


Then,  Interpolating 

'’c  ■ 2 - B)(Qa  + %)  + BQ, 

for  values  at  point  c, 


Wp  - (1  - e)(w^  + wj  + pwg 

"'c  “ 2 ^d) 
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Let 


K, 


■■A 

Introduce  nondlraenalonal  variables 


JTa 


JL 

pAcJ 


N„ 

pAc. 


M. 


At 

Define  the  following  constants: 
Y 


'^A 

% 

' % 

\h 

picj 

plCb 

A. 

pAc2 

n . A. 
B , 2 

pACg 

Wa 

“b 

w 

° % 

=b 

A. 

pAc^ 

pACb 

AtCb 

^A 

Al 

2 


1 pAK^ 


V - ^ 
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- -a3  - 4v 

bj  - + yj/b^  - (1  + 0)/[y(3  - g)] 

^3  - (1  + g)/[gv(3  - 0)] 

1^4  - ajbj  - 4e^b3 

'i  ■ ■ \ ■ 

■[<1  - 8)(!!„  + 17^/9  . + 29^5^  . p), 

^ The  .olutlon  to  I,.,  (l).(6)  to  t„„,  i. 

■ [»i(’'b  + ’b  • '^d)  • ‘vF2^]/(«2  + 4v) 

'^A  -"b 

''a  ■ (“■iCi  * - AvbjFj/b^ 

"a  ■ ("b'Ta  ■ cJ/kj 

\ ■ |-»3®D  + ‘vfH,  - 3^  - 9“4iJ„  - B^(.2  + v)5jj/b; 

®A  ■ - »3(®*  - 

SlmHarly,  the  equations  for  the  rear  end  of  the  beam  are 

"b  ■ "x  • “Vz^B  - "x)  ' 0 (! 

"Z  - “x  + ‘■Iz'bK  - •*)  - s(^^^]At  (8 
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pAc. 


- ”z 
^ % 


. -JL 

X . 2 

pACg 

Wm 

ty  . Ji 

% 


UJ- 


At 


z z 
% 

Ate. 


1 ^ 
^ 'b 


Define : 


V - ^ 

^ 2 

"a 

"4  pAKj 

Qlm  • ~ 

5 

®6  " plK^ 

■ l/[v(3  - p)] 

^2  “ dj^/p 
^3  - “ttg  - 

- l/[2(v+  Oj)] 
d;  ■ (2  - d3d,)/4 
dg  - -0^V/2  - (1  + a)dj^d5 

d^  - (1  + P)d2d3 
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+ “s'l 


+ 2pijJ  + d2[<i  - »)%  * 2pnJ 

+[(1  - + 2puij^]/(3  - 0) 

• '*4[»5K  *\*  + “s'l 

The  solution  to  Eqs,  (7) -(12)  in  terns  of  dimensionless  quantities  Is: 
" “4(\  ■ \)/(®4 


- i?2  " \ 


^Z  " (®2  ■ ^6^z)^^7 

uJ^  - (1  + 9)  {^2  ‘ ‘ ®1 

^Z  “ • ®6(\  " ®z)/(4y) 
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Laurence  Livermore  Laboratory 
ATTNi  Jerry  Coudreau 
ATTNi  Tech,  Info.,  Dept.  L-.1 
ATTNi  Mark  Wllklna,  L-504 


Loa  Alamoa  Scientific  Laboratory 

ATTNi  Doc.  Control  for  Reporta  Lib. 
ATTNi  Doc.  Control  for  Tom  Dowlor 


Livermore  Laboratory 
Control  for  Tech,  Library 


Sandla  Laboratoriea 
ATTNi  Doc. 

Sandia  Laboratoriea 


ATTNi 

Doc  s 

Con. 

ATTNi 

Doc. 

Con. 

ATTNi 

Doc . 

Con, 

ATTNi 

Doc . 

Con. 

ATTNi 

Doc . 

Con. 

ATTNi 

Doc . 

Con. 

ATTNi 

Doc . 

Con. 

For  John  Colp 

For  John  Keltur 

For  William  Patteraon 

for  William  Caudle 

For  W.  Altameirer 

for  3141  Sandla  Rpt.  Coll 

for  Walter  Herrmann 


EOiC,  Inc. 

Albuquerque  Diviaion 

ATTNi  Technical  Library 

Engineering  Socletlea  Library 
ATTNi  Ann  Mott 

Oeneral  Dynaraici  Corp. 

Pomonii  Division 

ATTNi  Keith  Anderson 

General  Electric  Company 

TEMPO-Center  for  Advanced  Studies 
ATTNi  DA51AC 

Cleorglo  Institute  of  Technology 

Oeorgln  Tech.  Research  Inatitutu 
ATTNi  L,  W.  RehField 
ATTNi  S,  V.  Hanagud 

Honeywell  Incorporated 

Defense  Syatema  Dlviai  on 
ATTNt  Ts  Ns  HbIvIr 
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Inicltut*  for  DtftnM  Analyiai 

ATTNl  IDA  Librarian,  Ruth  S.  Smith 

Kaman  AvlDyna 

Division  of  Katiian  Sclancaa  Corpt 
ATTNl  E.  S,  Crlieiona 
ATTNl  Norman  V,  Hobba 
ATTNl  Tachnioal  Library 

Kaman  Sclancaa  Corporation 
ATTNl  Library 

Lockhcad  Hlsallaa  & Spaca  Co.,  Inc. 

ATTNl  Tachnioal  Library 
ATTNl  M.  Culp 

Lockhaad  Hlasllaa  and  Spaca  Co.,  Inc. 

ATTNl  Taoh,  Info,  Ctr.,  D/Coll. 

Martin  Karlatta  Corporation 

Orlando  Dlvlalon 

ATTNl  A1  Cowan 
ATTNl  M,  Anthony 
ATTNl  H.  McQuatg 

Karrltt  CASES,  Incorporatad 
ATTNl  J,  L.  Marrltt 
ATTNl  Technical  Library 

Untvaralty  of  Naw  Haaloo 

Oapt.  of  Campus  Security  and  Police 
ATTNl  a.  E.  Trlandafalldla 

Nathan  M.  Nawmark 

Consulting  Engineering  Services 
ATTNl  W.  Hall 
ATTNl  Nathan  M.  Nawmark 

Pacifica  Taohnology 
ATTNl  R,  BJork 
ATTNl  0,  Kant 

Physics  Intarnatlonnl  Company 

ATTNl  Doc.  Con.  for  Tech,  Lib. 

ATTNl  Doc.  Con.  for  Dennis  Orphnl 
ATTNl  Doc.  Con.  for  Larry  A Hehrmnnn 


DEPARTMENT  OF  DEFENSE  CONTRACTOns 


R t 


D Aasoclatas 


ATTNl 

ATTNl 

ATTNl 

ATTNl 

ATTNl 

ATTNl 

ATTNl 

ATTNl 


Cyrus  P.  Knowlas 
J.  0.  Lewis 

William  B,  Wright,  Jr. 
Paul  Rauaoh 
Harold  L.  Brode 
Tachnlcal  Library 
Arlan  Plaids 
Hanry  Coops r 


Tha  Rand  Corporation 

ATTNl  Tachnlcal  Library 


Sclanca  Applications,  Inc. 

ATTNl  Tachnlcal  Library 


SRI  International 

ATTNl  Oeorge  R,  Abrahamaon 
ATTNl  Jim  Colton 
ATTNl  H.  E.  Llndberg 
ATTNl  R.  E.  Emeraon 
ATTNl  J.  K.  Oran 


Syatama,  Sclanca  and  Software,  Inc. 
ATTNl  Edward  Oaffnay 
ATTNl  Hobart  Sadgewlck 
ATTNl  Technical  Library 

Terra  Tak,  Inc, 

ATTNl  Tachnlcal  Library 

TRW  Defenaa  & Space  Sya.  Croup 

ATTNl  Tech.  Info.  Oanter/S-lSlO 
ATTNl  Pater  K.  Dal,  Rl/2170 


TRW  Defense  i.  Space  Bys,  Croup 
San  Bernardino  Operatlona 


ATTNl 

E.  y.  Wong,  527/712 

Weidllnger  Aaaoc.  Conaultlng  Enalnecra 

ATTNj 

M.  McCormick 

ATTNt 

Melvin  L,  Boron 

t/cidUnger  t 

^iioc.  ConiUltlnK  Etiglneora 

ATTNt 

vK  tgenberg 

I 

I 
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